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^  Effective,  comprehensive,  and  low  cost  procedures  do  not  exist  for 
eliminating  scour  during  construction  in  the  nearshore  environment.  Deter¬ 
mination  of  potential  alternative  procedures  Is  seriously  hampered  by  the 
inability  to  predict  the  extent  of  potential  scour.  The  objectives  of  this 
research  program  are  to  develop  techniques  to  minimize  and  control  scour 
during  nearshore  construction  and  to  predict  the  probable  magnitude  of 
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20.  ABSTRACT  (Continued) 

r-"  scour  that  may  result  as  a  function  of  currents  and  wave  climate.  One 
phase  of  the  research  effort  is  the  development  of  numerical  techniques 
(incorporating  both  refraction  and  diffraction  effects  near  the  structure) 
for  computing  wave-induced, velocities  and  tidal  currents  in  the  vicinity  . 
of  structures  and  applying  these  results  to  determine  sediment  transport 
of  the  bottom  material  at  the  particular  site. 

The  purpose  of  this  study  was  to  obtain  detailed  and  precise 
experimental  data  regarding  wave-height  variations  and  currents  (patterns 
and  magnitudes)  downwave  from  a  shore-connected  breakwater  or  jetty  under 
the  simultaneous  effects  of  refraction  and  diffraction.  This  information 
provides  insight  into  the  phenomenon  of  combined  wave  refraction  and 
diffraction  and  can  be  used  to  verify  numerical  models  that  simulate  this 
phenomenon^  j  The  experimental  investigation  was  conducted  in  a  wave  basin 
that  was  molded  in  cement  mortar  and  consisted  of  an  area  50  ft  x  60  ft 
with  a  water  d^pth  of  1  ft  in  the  open-ocean  region.  A  vertical,  im¬ 
permeable  breakwater  (shore-connected)  extending  perpendicularly  from  the 
shoreline  was  installed  on  a  beach  slope  of  IV  on  20H. 

In  the  neighborhood  of  the  breakwater,  currents  existed  that  af¬ 
fected  the  wave  heights.  The  magnitude  of  these  wave-induced  return  cur¬ 
rents  is  a  function  of  incident  wave  characteristics.  The  effect  of  varying 
the  incident  wave  height  on  the  wave-height  amplification  factor,  H/H  , 
was  investigated  and  it  was  determined  that  the  greatest  variation  in  H/H 
occurred  in  the  deep  shadow  zone  near  the  breakwater  and  shore  region 
where  the  currents  are  the  strongest.  This  nonlinear  effect  diminishes 
rapidly  away  from  the  structure,  and  at  a  distance  of  5  ft  (model  dimen¬ 
sions)  downwave  appears  to  be  relatively  insignificant  except  for  the  ex¬ 
treme  shadow  region  at  acute  angles  of  incidence.  The  establishment  of 
a  counterclockwise  circulation  cell  approximately  4  ft  wide  adjacent  to 
the  downwave  side  of  the  breakwater  results  in  a  seaward  flowing  current 
all  along  the  breakwater.  The  bottom  current  is  especially  intense 
approximately  40  percent  of  the  breakwater  length  from  the  shoreline  and 
decreases  seaward  along  the  structure.  The  physical  model  sidewall 
boundary  was  responsible  for  a  clockwise  circulation  cell  which  developed 
as  a  result  of  mass  transport  downcoast  by  the  nearshore  current.  Based 
on  these  experiments,  any  numerical  techniques  for  describing  scour  and 
erosion  near  structures  should  account  for  wave-induced  currents  and 
circulation  cells  near  the  structure  as  well  as  tidal  or  other  currents 
that  may  exist  in  the  near  region. 
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PREFACE 


The  study  reported  herein  was  authorized  as  a  part  of  the  Civil 
Works  Research  and  Development  Program  by  the  Office,  Chief  of  Engineers 
(OCE).  This  particular  work  unit.  Erosion  Control  of  Scour  During  Con¬ 
struction,  is  part  of  the  Improvement  of  Operations  and  Maintenance 
Techniques  (IOMT)  Program*  Mr.  William  Godwin  was  the  OCE  Technical 
Monitor  for  the  IOMT  Program  during  preparation  and  publication  of  this 
report. 

This  study  was  conducted  during  the  period  1  March  1979  through 
30  April  1980  by  personnel  of  the  hydraulics  Laboratory  of  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  under  the  general  supervision 
of  Messrs.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory;  F.  A. 
Herrmann,  Jr.,  Assistant  Chief  of  the  Hydraulics  Laboratory;  R.  A.  Sager, 
Chief  of  the  Estuaries  Division  and  IOMT  Program  Coordinator;  Dr.  R.  W. 
Whalin,  Chief  of  the  Wave  Dynamics  Division;  Mr.  D.  D.  Davidson,  Chief 
of  the  Wave  Research  Branch;  and  Dr.  J.  R.  Houston,  Research  Engineer 
and  Principal  Investigator  for  the  Erosion  Control  of  Scour  During  Con¬ 
struction  work  unit.  Dr.  L.  Z.  Hales,  Research  Hydraulic  Engineer, 

Mr.  K.  A.  Turner,  Computer  Specialist,  Ms.  M.  L.  Hampton,  Computer  Tech¬ 
nician,  Mr.  R.  E.  Ankeny,  Computer  Technician,  and  Mr.  K.  M.  Strausbaugh, 
Civil  Engineering  Technician,  performed  the  study  described  herein. 

Dr.  Hales  prepared  this  report. 

Commanders  and  Directors  of  WES  during  the  conduct  of  this  in¬ 
vestigation  and  the  preparation  and  publication  of  this  report  were 
COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE.  Technical 
Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (SI)  units  as  follows: 

Multiply _  _ By  _  To  Obtain 


feet  0. 

feet  per  second  0. 

feet  per  second  per  0* 

second 

foot-pounds  per  foot  148. 

per  foot 

foot-pounds  per  second  13825. 

Inches  25. 

pounds  per  foot  per  4, 

foot 

pounds-second  per  foot  4. 

per  foot 

pounds-second-second  per  52. 

foot  per  foot  per  foot 
per  foot 

square  feet  0. 


3048 

metres 

3048 

metres  per  second 

3048 

metres  per  second  per 
second 

8164434 

kilogram-centimetre 
per  metre  per  metre 

5 

gram-centimetre  per 
second 

4 

millimetres 

8824276 

kilogram  per  metre  per 
metre 

8824276 

kilogram-second  per 
metre  per  metre 

5540137 

ki log  ram-second-second 
per  metre  per  metre 
per  metre  per  metre 

09290304 

square  metres 

EROSION  CONTROL  OF  SCOUR  DURING  CONSTRUCTION 


EFFECT  OF  SIMULTANEOUS  REFRACTION  AND  DIFFRACTION 
ON  WAVE  HEIGHTS  AND  CURRENT  PATTERNS 


PART  I:  INTRODUCTION 


Statement  of  the  Problem 


1*  When  major  structures  are  erected  in  the  coastal  zone,  they 
alter  currents  that  are  in  a  dynamic  equilibrium  with  the  existing 
bathymetry.  These  altered  currents  may  change  the  existing  bathymetry. 

In  addition,  waves  breaking  on  the  new  structure  will  cause  bottom  ma¬ 
terial  to  be  suspended  and  transported  from  the  region  by  longshore  or 
other  currents.  This  removal  of  material  from  around  structures  is 
often  not  compensated  by  an  influx  of  additional  material;  the  result  is 
scour,  or  erosion,  that  usually  develops  along  the  toe  of  the  structure. 
In  order  to  ensure  structural  stability,  the  scour  area  must  be  filled 
with  nonerodible  material  (sufficiently  stable  to  withstand  the  environ¬ 
mental  forces  to  which  it  will  be  subjected).  This  may  result  in  addi¬ 
tional  quantities  of  material  being  required  during  construction  that 
can  potentially  be  very  costly.  To  minimize  potential  cost  increases 
due  to  scour  during  construction,  it  is  necessary  to  quantify  the 
probability  and  ultimate  extent  of  potential  scour  during  the  scheduled 
construction  period. 

2.  Effective,  comprehensive,  and  low  cost  procedures  do  not  exist 
for  eliminating  scour  during  construction  in  the  nearshore  environment. 
Determination  of  potential  alternative  procedures  is  seriously  hampered 
by  the  inability  to  predict  the  extent  of  potential  scour.  Objectives 
of  this  research  effort  are  to  develop  techniques  to  minimize  and  control 
scour  during  nearshore  construction  and  to  predict  the  probable  magnitude 
of  scour  that  may  result  as  a  function  of  currents  and  wave  climate. 

One  phase  of  the  research  program  is  development  of  numerical  techniques 
(incorporating  both  refraction  and  diffraction  effects  near  the 
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structure)  for  computing  wave-induced  velocities  and  tidal  currents  in 
the  vicinity  of  structures  and  applying  these  results  to  determine 
sediment  transport  of  the  bottom  material  at  the  particular  site.  This 
will  result  in  an  accounting  of  the  flux  of  material  around  the  struc¬ 
ture,  and  thus  knowledge  of  the  extent  of  erosion  or  accretion  to  be 
expected  as  a  function  of  wave  climate,  currents,  local  bathymetry,  and 
structure  design. 

Purpose  of  the  Study 

3.  The  present  state  of  nearshore  current  and  wave  theories  has 
reached  the  point  where  detailed  experimental  investigations  are  re¬ 
quired  for  the  verification  of  analytical  developments  and  numerical 
models.  To  provide  a  firm  foundation  for  further  advancements,  a  simple 
beach  profile  consisting  of  straight,  uniform  contours  parallel  with 
the  shoreline  was  physically  modeled.  A  shore-connected  vertical,  im¬ 
permeable  barrier  (breakwater)  was  installed  perpendicular  to  the  shore¬ 
line  to  simulate  the  usual  prototype  jetties  and  breakwaters  commonly 
occurring  along  many  coasts.  This  single  jetty  (shore-connected  break¬ 
water)  simplifies  the  experiment,  facilitates  direct  comparisons  with 
numerical  model  results,  and  provides  greater  understanding  and 
insight  into  the  phenomena  of  wave-height  variations  and  current  cir¬ 
culation  patterns  downwave  of  a  jetty  or  breakwater.  Analytical  models 
of  this  physical  arrangement  have  been  developed  (for  example,  Liu  and 
Mei  1975)  that  extrapolate  the  small  amplitude  wave  theory  to  calculate 
the  radiation  stress  distribution.  Such  analytical  developments  re¬ 
garding  practical  cases  involving  both  refraction  and  diffraction  have 
not  been  verified  for  either  wave-height  distribution  or  currents 
(magnitude  and  direction)  in  the  near-vicinity  of  the  structure. 
Knowledge  of  these  important  phenomena  will  be  used  as  the  basis  for  ad¬ 
vanced  studies  of  sediment  movement  around  major  structures  under  com¬ 
bined  effects  of  wave  refraction  and  diffraction. 

4.  The  purpose  of  this  study  was  to  obtain  quantitative  and  de¬ 
tailed  laboratory  measurements  of  combined  refraction  and  diffraction  in 


the  lee  of  a  jetty  or  shore-connected  breakwater.  These  measurements 
include  both  the  wave-height  distribution  downwave  of  the  structure,  and 
the  associated  wave-induced  current  field  and  circulation  patterns,  and 
can  be  used  to  verify  numerical  models  that  simulate  the  phenomena. 
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PART  II:  LITERATURE  REVIEW 


5.  The  approximate  solution  of  water  wave  refraction  caused  by  a 
variable  bathymetry  Is  well  known  and  can  be  derived  by  assuming  that 
bottom  reflections  are  negligible  or  by  the  more  rigorous  Wentzel- 
Kramers-Brillouin  (WKB)  approximation.  The  exact  solution  for  the  dif¬ 
fraction  of  surface  waves  by  vertical  barriers  of  simple  cross  section 
in  water  of  constant  depth  also  is  well  known  (being  analogous  to  clas¬ 
sical  problems  of  physics).  An  analytical  theory  for  the  practical  case 
of  combined  refraction  and  diffraction  has  not  been  completely  developed. 
The  present  engineering  practice  for  determining  wave  heights  under  this 
condition  is  a  stepwise  procedure  (Dunham  1951;  Liu  and  Mei  1976a;  U.  S. 
Army  Engineer  Coastal  Engineering  Research  Center*  CERC  1977).  The  pro¬ 
cedure  involves  the  following  steps:  (a)  calculate  refraction  effects  up 
to  the  barrier  (jetty  or  shore-connected  breakwater) *  (b)  calculate  for  a 
"few"  wavelengths  the  effects  of  diffraction  assuming  a  constant  water 
depth*  and  (c)  beyond  this  region  calculate  the  refraction  effects  only. 
This  procedure  is  obviously  imprecise,  and  Mobarek  (1962)  indicates  that 
the  method  is  suitable  only  for  intermediate  water  depths.  In  addition* 
it  can  be  stated  that  this  procedure  is  only  valid  for  small  refraction 
effects  as  demonstrated  by  Whalin  (1972). 

6.  In  addition  to  wave-height  variation,  scour  and  erosion  around 
structures  are  also  influenced  significantly  by  current  magnitudes  and 
circulation  patterns.  These  phenomena  can  develop  differently  under  the 
simultaneous  effects  of  refraction  and  diffraction.  Many  studies* 
theoretical  and  laboratory*  have  investigated  the  uniform  longshore 
current  on  a  long  beach  of  straight,  parallel  contours,  since  this  par¬ 
ticular  current  is  many  times  the  roost  important  of  littoral  transport 
mechanisms.  Early  analytical  work  considered  the  longshore  component 

of  the  momentum  flux  to  be  the  net  driving  force  to  be  balanced  by 
bottom  friction  (proportional  to  the  square  of  the  longshore  velocity). 
This  resulted  in  an  indication  of  the  magnitude  of  the  longshore  current. 
The  advances  of  Longuet-Higgins  and  Stewart  (1960*  1961,  1964),  Longuet- 
Higgins  (1970a,  1970b),  and  Lundgren  (1962)  were  important  because  of 
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the  introduction  of  the  radiation  stress.  This  stress  is  induced  by  wave 
fluctuations  and  is  the  time  average  of  the  local  horizontal  momentum 
flux  integrated  over  the  depth.  Away  from  the  surf  zone,  these  stresses 
have  a  sound  theoretical  basis.  In  the  surf  zone,  where  turbulence  due 
to  breaking  conditions  exists,  theories  for  describing  either  the  oscil¬ 
latory  wave  motion  or  radiation  stress  have  not  been  fully  developed. 
Hence,  the  relationships  between  radiation  stress  and  mean  flow  veloc¬ 
ities  have  been  hypothesized  only  by  empiricism  and  plausible  reasoning. 

7.  Theories  have  been  developed  by  Longuet-Higgins  (1970a,  1970b) 
and  Thornton  (1970)  for  describing  longshore  currents  on  a  straight  beach 
with  parallel  contours  which  rely  on  the  radiation  stress  concept.  Em¬ 
pirical  reasoning  regarding  lateral  turbulence  and  bottom  friction 
allows  reasonable  velocity  distributions  to  be  obtained  for  the  longshore 
current.  Liu  and  Mei  (1975,  1976b)  have  investigated  analytically  the 
combined  effects  of  refraction  and  diffraction  on  both  wave  heights  and 
currents  due  to  a  jetty  or  shore-connected  breakwater  on  a  sloping  beach 
by  ignoring  convective  inertia  and  lateral  turbulent  diffusion.  The 
resulting  equations  are  solved  by  finite  difference  techniques  and,  in 
addition  to  the  displacement  of  the  longshore  current,  a  cell  develops 
in  the  shadow  zone.  Vorticity  arguments  indicate  that  the  radiation 
stress  establishes  a  mean  sea  level  gradient  which  drives  the  longshore 
current • 


Refraction  and  Shoaling 

8.  In  intermediate  and  shallow  water,  the  phase  speed  of  a  surface 
gravity  wave  depends  on  water  depth.  Since  wave  celerity  decreases  as 
depth  decreases,  the  wavelength  must  also  decrease  for  the  period  to 
remain  constant.  Phase  velocity  varies  along  the  crest  of  a  wave  propa¬ 
gating  at  an  angle  to  underwater  contours  because  that  part  of  the  wave 
in  deeper  water  moves  faster  than  that  part  in  shallower  water.  This 
variation  causes  the  wave  crest  to  bend  toward  alignment  with  the  con¬ 
tours.  This  bending  effect  due  to  changes  in  bottom  topography,  called 
refraction,  depends  on  the  relation  of  water  depth  to  wavelength  and  is 


8 


analogous  to  refraction  of  other  types  of  waves  such  as  light.  A  basic 
assumption  in  wave  refraction  theory  is  the  conservation  of  energy 
between  wave  orthogonals  (i.e.,  no  diffraction  of  energy  along  wave 
crests) . 

9.  As  waves  propagate  from  deep  water  into  shallow  water,  changes 
other  than  refraction  occur.  These  changes,  called  shoaling,  are  in¬ 
vestigated  under  the  assumption  that  there  is  no  loss  of  wave  energy  and 
negligible  reflection.  The  power  being  transmitted  by  the  wave  train  in 
water  of  any  depth  is  equal  to  the  power  being  transmitted  by  the  wave 
system  in  deep  water.  The  wave  period  remains  constant  regardless  of  the 
water  depth,  whereas  the  wavelength,  velocity,  and  height  vary. 

10.  Refraction  and  shoaling  effects  are  important  for  several 
reasons.  These  phenomena  determine  the  wave  height  in  many  locations  for 
a  given  set  of  incident  deepwater  wave  conditions;  i.e.,  wave  height, 
period,  and  direction  of  propagation  in  deep  water.  Refraction  and 
shoaling,  therefore,  have  a  significant  influence  on  the  distribution 

of  wave  energy  along  the  coast.  The  change  in  wave  direction  of  differ¬ 
ent  parts  of  the  wave  results  in  convergence  or  divergence  of  wave  energy 
and  materially  affects  the  forces  exerted  by  waves  on  structures  and  of 
the  capacity  of  waves  to  transport  sand  either  alongshore  or  onshore/ 
offshore. 

Waves  on  water  of  variable  depth 

11.  Procedures  for  the  computation  of  refraction  of  surface 
gravity  waves  on  water  of  nonuniforra  depth  involve  the  assumption  that 
a  wave  with  a  curved  crest  pattern  and  variable  amplitude  along  the 
crest  behaves  locally  as  a  straight-crested  wave  of  constant  amplitude. 
Rayleigh  (1877)  appears  to  have  been  the  first  to  use  the  approximations 
of  geometrical  optics  in  this  analysis,  and  theoretical  results  have 
been  developed  with  respect  to  energy  flux  and  phase  speed.  As  expressed 
by  Keller  (1958),  the  geometrical  optics  theory  defines  a  propagation 
velocity  at  each  point  on  the  water  surface,  with  this  velocity  being 
exactly  that  which  waves  of  given  period  would  have  in  water  of  uniform 
depth  at  all  points.  By  employing  Fermat’s  principle  of  optics,  wave 
rays  are  defined  and  surface  waves  are  assumed  to  propagate  along  these 


rays.  The  variation  of  the  amplitude  along  the  rays  is  determined  by 
the  use  of  the  principle  of  conservation  of  energy*  This  principle  (in 
its  optical  form)  states  that  the  flux  of  energy  is  the  same  at  all  cross 
sections  between  two  adjacent  wave  rays*  The  energy  flux  is  proportional 
to  the  square  of  the  amplitude  of  the  waves  and  to  the  distance  between 
the  rays*  and  hence  the  wave-height  variation  along  the  ray  is  available, 

12,  Keller's  (1958)  derivation  is  based  upon  an  asymptotic  solu¬ 
tion  of  the  equations  of  the  exact  linear  theory  for  periodic  waves  in 
water  of  arbitrary  nonuniform  depth.  The  solution  is  asymptotic  in  the 
sense  that  the  depth  and  wavelength  are  small  compared  with  the  hori¬ 
zontal  scale  of  the  bottom  contours.  The  first  term  of  the  solution 
agrees  exactly  with  the  asymptotic  form  of  the  solution  for  waves  in 
water  with  a  uniformly  sloping  bottom  as  the  bottom  slope  tends  to  zero* 
This  solution  conforms  with  all  the  principles  of  the  geometrical  optics 
theory  of  Rayleigh  (1877)  and  thus  provides  a  derivation  of  that  theory* 
The  results  are  not  valid*  however*  at  caustics  or  ray  crossings, 

13.  In  problems  of  linear  wave  propagation  over  mild  slopes,  the 
principle  of  geometrical  optics  has  been  applied  by  Carrier  (1966)  as 
the  first  approximation  in  a  systematic  perturbation  scheme  while  the 
bottom  is  considered  to  be  locally  horizontal.  The  depth  variation  was 
dealt  with  afterward  by  requiring  the  appropriate  energy  conservation. 
This  was  also  the  basic  idea  for  the  work  of  Koh  and  Le  Mehaute  (1966) 

in  which  the  transformation  of  progressive  waves  was  investigated  as  they 
travel  from  deep  water  to  shore.  The  Stokes'  theory  at  a  fifth  order  of 
approximation  was  applied  along  with  the  method  of  conservation  of 
energy  flux.  The  first,  third*  and  fifth  orders  of  approximation  were 
compared  with  each  other  and  with  experiments*  The  differences  between 
the  predictions  of  wave-height  changes  based  on  the  three  orders  of 
approximation  were  found  to  be  small*  on  the  order  of  5  percent.  For 
practical  purposes,  the  third-order  theory  was  found  to  give  reliable 
results*  The  third  and  fifth  order  Stokes’  theories  are  based  on  a 

3 

series  expansion  in  terms  of  H/L  where  terms  of  the  order  of  (H/L) 
and  (H/L)**  ,  respectively,  are  retained  and  higher  order  terms  are 
neglected.  It  should  be  noted  that  this  theory  is  based  on  an  expansion 
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* 


in  term  of  the  wave  steepness,  H/L  ,  and  consequently  can  be  expected 
to  better  approximate  limit  steepness  waves  in  deep  water.  However,  it 
cannot  be  expected  to  do  a  very  good  job  of  approximating  waves  in  shal¬ 
low  water  since  water  depth  is  not  a  parameter  in  the  series  expansion. 

It  also  is  assumed  that  the  wave  is  simply  harmonic  in  time.  This 
theory  could  be  considered  as  a  finite  amplitude  deepwater  wave  theory. 

14.  In  cases  of  limiting  shallow  water,  the  wave  conditions  are 
nearly  independent  of  wavelength,  and  the  important  parameters  are  water 
depth  and  the  ratio  of  wave  height  to  water  depth.  Mathematical  argu¬ 
ments  show  that  Stokes  waves  are  most  nearly  valid  in  water  deeper  than 
about  d/L  >  1/8  to  1/10  (Keulegan  1950).  In  shallower  water,  the  the¬ 
ory  for  a  wave  type  known  as  cnoidal  appears  to  be  more  satisfactory, 
and  Masch  (1964)  Investigated  the  problem  of  wave  shoaling  using  cnoidal 
wave  theory  with  the  formulas  developed  by  Keulegan  and  Patterson  (1940). 
Masch  (1964)  assumed  hydrostatic  pressure  distribution  and  neglected  the 
convective  inertia  term  in  his  expression  for  the  energy  flux.  The 
third  and  fifth  order  cnoidal  theories  are  based  on  a  series  expansion 

in  terms  of  H/d  ,  where  terms  of  the  order  of  (H/d)  and  (H/d)  , 

respectively,  are  retained  and  higher  order  terms  neglected.  It  should 
be  noted  that  this  theory  is  based  on  an  expansion  of  the  relative  wave 
height  (H/d)  and  can  be  expected  to  better  approximate  the  wave  form  in 
shallow  water.  However,  it  cannot  be  expected  to  do  a  very  good  job  of 
approximating  the  wave  form  for  limiting  steepness  waves  in  deep  water. 

In  that  case,  water  depth  is  unimportant  and  wavelength  Is  crucial.  This 
theory  could  be  considered  a  finite  amplitude  intermediate  and  shallow- 
water  theory. 

15.  A  technique  of  asymptotic  expansion  was  applied  by  Mei,  Tlapa, 
and  Eagleson  (1968)  to  water  wave  propagation  over  an  uneven  bottom  that 
has  straight  and  parallel  contours.  Attention  was  focused  on  the  estab¬ 
lishment  of  a  rigorous  scheme  of  successive  approximation  for  higher 
order  corrections.  The  bottom  depth  was  assumed  to  vary  slowly  within  a 
wavelength.  By  introducing  a  compressed  coordinate  in  the  direction 
normal  to  the  contours  and  by  assuming  an  expansion  of  the  WKB  type,  the 
weakness  of  the  depth  variation  in  the  normal  direction  was  incorporated 
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In  the  mathematical  formulation.  The  conventional  linearized  theory  of 
wave  refraction  was  obtained  as  the  first-order  solution  without  the 
explicit  assumption  of  Snell's  law.  In  the  second  order,  a  steady-state 
depression  of  the  mean  water  surface  was  found  for  the  general  case  where 
the  incident  wave  approaches  the  contours  obliquely.  To  the  first  order, 
this  development  was  the  same  as  the  classical  theory  of  Rayleigh  (1877). 
At  higher  order  some  differences  existed  from  the  usual  adaptation  of  the 
Stokes'  theory  for  a  horizontal  bottom.  The  source  of  the  difference  is 
the  explicit  appearance  of  the  variation  of  the  bottom  boundary  condi¬ 
tions  for  the  first  harmonic  at  second  and  higher  orders.  The  theory  is 
expected  to  hold  up  to  the  neighborhood  where  the  wave  breaks,  except  at 
the  shoreline  where  a  singularity  exists  for  all  orders. 

16.  The  refraction  of  surface  gravity  waves  propagating  in  an 
ideal  fluid  was  investigated  by  Battjes  (1968)  with  amplitudes  small  so 
that  linear  theory  was  applicable  and  harmonic  in  time.  It  was  known 
a  priori  that  the  velocity  of  propagation  of  a  wave  crest  (to  third  order 
and  greater)  is  a  function  of  the  wave  height.  In  zones  of  convergence 
or  divergence  of  wave  energy,  gradients  in  wave  height  will  exist  along 
sections  of  the  wave  crest.  The  regions  of  greater  wave  height  will 
propagate  faster  than  the  regions  of  lesser  wave  height  and  this  will, 
in  turn,  create  bending  of  the  wave  crest  (in  addition  to  that  bending 
caused  by  the  bottom  topography,  called  refraction).  This  supplemental 
bending  is  not  usually  accounted  for  by  refraction  analysis.  Battjes 
(1968)  developed  an  expression  for  this  local  correction  to  the  wave 
speed  because  of  wave-height  variations  along  the  crest  from  an  exact 
derivation  of  the  wave  number,  k  .  It  was  determined  that  the  magnitude 
of  the  exact  wave  number  equals  k  plus  corrections  which  depend  on  the 
second  derivative  of  the  wave  amplitude.  This  correction  of  the  wave 
number  implies,  in  turn,  a  correction  to  the  phase  speed  c  ,  which  was 
found  to  be: 


c 


(1) 
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where 


2 

g  -  gravitational  constant,  32.174  ft/sec  * 
k  -  wave  number,  2tt/L,  i/ft 
d  •  local  water  depth,  ft 
a  ■  local  wave  amplitude,  ft 

The  second  derivative  of  the  wave  amplitude  in  the  horizontal  plane  is 

given  by  a  or  a  The  rate  of  power  transmission  P  ,  or  energy 

xx  yy 

flux,  was  determined  to  be: 


PAb  -  (1/2  pga2n  ~  Ab)|l  ♦  XX 


/  »  “  V'2 

l1*^) 


(2) 


where 


Ab  *  wave  ray  spacing,  ft 

2  4 

p  -  fluid  density,  lb-sec  /ft 

n  ■  ratio  of  group  velocity  c  ,  to  phase  velocity  c  , 
dimensionless  ® 

^  *  angular  frequency,  2tt/T  ,  1/sec 


Define: 


(3) 


The  commonly  used  existing  procedures  (Dunham  1951,  Liu  and  Mei  1976a, 
CERC  1977),  for  construction  or  computation  of  refraction  diagrams 
utilize  phase  speeds  that  are  obtained  by  neglecting  6  . 

17.  Battjes  (1968)  examined  the  omission  of  6  from  a  qualitative 
standpoint.  In  an  area  of  strong  local  convergence,  omission  of  6 
from  Equation  1  results  in  underestimating  the  local  phase  speed.  The 
result  is  that  in  Equation  2  the  ray  separation  Ab  will  be  under¬ 
estimated.  Thus,  there  are  generally  two  contributions  to  the  error 
which  results  in  the  computed  wave  amplitude  a  •  However,  these  two 
contributions  are  of  a  different  nature  because  the  effect  on  the  wave 
pattern  of  using  an  incorrect  value  for  the  phase  speed  is  cumulative. 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ments  to  metric  (SI)  units  is  presented  on  page  3. 


k 


whereas  the  effect  on  energy  flux  is  local.  In  any  case,  the  omission 
of  6  will  generate  wave  heights  at  variance  with  the  height  inferred 
from  refraction  diagrams  based  on  linear  small  amplitude  wave  theory 
which  neglects  the  effects  of  wave-height  gradients  along  the  wave 
crest,  so  that  it  would  appear  that  energy  had  been  transferred  across 
orthogonals.  An  estimate  of  the  magnitude  of  8  has  been  approximated 
by  Battjes  (1968)  for  four  distinctly  different  cases.  For  simple  shoal¬ 
ing,  the  amplitude  varies  only  in  the  direction  of  wave  propagation,  and 
for  shallow  water  the  variation  was: 


ft  -  8  x  10-3  (L/d)2S2 


where  S  is  the  bottom  slope.  For  short-crested  waves  where  the 
distance  along  the  wave  crest  is  two  or  three  times  the  wavelength, 

8  -  -20%  or  -10%  ,  respectively.  For  the  case  of  diffraction  around  a 
semi-infinite  breakwater,  8  was  found  to  reach  values  between  +10% 
and  -10^  at  points  a  distance  of  one  wavelength  from  the  breakwater  tip, 
decreasing  Inversely  proportional  to  the  distance  from  the  tip.  For  the 
case  of  diffraction  through  a  gap  of  width  two  times  the  wavelength,  6 
was  evaluated  in  a  few  points  on  the  center  line  of  the  gap  where  it  was 
found  to  reach  values  of  25%,  7%,  and  3%  at  distances  of  1,  2,  and  4 
wavelengths,  respectively,  from  the  gap. 

Longshore  currents 
due  to  wave  breaking 

18.  It  is  well  known  that  as  sea  or  swell  approaches  a  coastline 
at  an  angle,  breaking  waves  create  a  mean  current  in  the  surf  zone  paral¬ 
lel  to  the  coastline.  This  longshore  current  transports  sand  or  other 
material  and  is  of  primary  Importance  to  coastal  engineers.  It  has 
often  been  suggested  that  the  magnitude  of  the  longshore  current  is  re¬ 
lated  in  some  way  to  the  energy  flux  or  momentum  flux  of  the  incoming 
wave  train.  Indeed,  it  appears  reasonable  to  believe  that  the  forces 
available  for  transporting  material  should  be  related  to  the  available 
wave  power  in  the  vicinity.  The  approach  employing  momentum  flux  should 
be  the  more  promising  since  momentum  is  conserved,  whereas  energy  is 
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dissipated  in  the  breaking  process.  Even  so,  Inman  and  Quinn  (1952) 
have  shown  that  in  order  to  make  the  momentum  flux  theory  fit  prototype 
observations,  the  friction  coefficient,  C  ,  would  have  to  vary  over  a 
range  of  3  to  4  orders  of  magnitude.  The  momentum  flux  theory  which  was 
analyzed  by  Inman  and  Quinn  (1952)  was  that  which  had  been  developed  by 
Putnam,  Munk,  and  Traylor  (1949)  and  contains  the  beach  slope  parameter. 

19.  Galvin  and  Eagleson  (1965)  performed  experimental  studies  of 
longshore  currents  on  a  plane  beach.  Measurements  were  made  of  the 
characteristics  of  breaking  waves  and  the  resulting  longshore  currents 
for  34  combinations  of  wave  height  (up  to  0.22  ft),  period  (0.90  to 
1.50  sec),  and  breaker  angle  (up  to  32  deg),  along  a  20-ft  test  section 
of  a  30-ft  plane,  smooth  concrete  beach  with  a  slope  of  0.104.  Observa¬ 
tions  and  measurements  showed  that  most  of  the  fluid  in  the  surf  zone 
stays  there,  and  that  longshore  current  velocity  Initially  increases 
downstream  from  an  obstruction  (jetty  or  shore-connected  breakwater). 
Velocity  Increases  along  the  beach  because  the  fluid  forming  the  breaking 
wave  has  been  withdrawn  from  the  surf  zone  and  thus  already  has  a  long¬ 
shore  component  of  motion  of  the  breaking  wave.  A  differential  equation 
for  this  nonuniform  flow  agreed  qualitatively  with  the  measured  variation 
of  velocity  with  breaker  angle  and  with  distance  from  an  obstacle.  The 
nonuniformity  of  the  flow  was  also  indicated  by  the  mean  water  level 
which  Increased,  and  by  the  breaker  position  and  runup  limit  which  moved 
shoreward  downstream  from  the  obstruction.  The  energy  used  to  maintain 
the  flow  of  the  longshore  current  was  found  to  be  a  small  fraction  (less 
than  10  percent  of  the  energy  brought  to  the  surf  zone  by  the  shoaling 
waves.  The  mean  velocity  of  the  uniform  longshore  current,  V^g  ,  was 
approximately: 


Vj  *  gST  sin  2 


(5) 


where 

S  -  beach  slope,  dimensionless 
T  •  wave  period,  sec 

0L  *  angle  breaking  wave  makes  with  the  shoreline,  deg 


Equation  5  Is  one  form  of  the  conservation  of  mass  in  the  surf  zone. 

20.  By  using  known  results  of  the  radiation  stress  associated 
with  gravity  waves.  Longue t-Higgins  (1970at  1970b)  investigated  the 
mechanism  of  obliquely  incident  sea  waves  which  generate  longshore  cur¬ 
rents.  The  total  lateral  thrust,  F^  ,  per  unit  distance  parallel  to 
the  coastline  exerted  by  incoming  waves  on  the  beach  and  in  the  nearshore 
zone  was  shown  to  be: 


F  -  1/4  E  sin  2  a 
W  o  o 


(6) 


where 

2 

Eq  ■  wave  energy  density  in  deep  water,  ft-lb/ft 
aQ  ■  deepwater  wave  approach  angle,  deg 
The  local  stress,  T  ,  per  unit  area  exerted  on  the  surf  zone  in  steady 
conditions  was  shown  to  be: 


x  *  D/c  sin  a 


(7) 


where 

D  -  local  rate  of  energy  dissipation,  lb/ft/sec 
c  -  phase  velocity,  ft/sec 
a  -  local  angle  of  incidence,  deg 

These  relations  appear  independent  of  the  manner  of  the  energy  dissipa¬ 
tion,  but  because  breaker  height  is  related  to  local  depth  in  shallow 
water,  it  is  conceivable  that  most  of  the  dissipation  is  due  to  wave 
breaking  and  not  to  bottom  friction.  Under  these  considerations,  the 
local  mean  longshore  stress  in  the  surf  zone,  x^g  ,  is  given  as: 

T,  -  5/4  p  uf  S  sin  a  (8) 

Is  max 

Here  u  is  the  maximum  orbital  velocity  in  the  wave,  and  the  other 
max 

symbols  have  been  previously  defined.  It  was  also  shown  that  if  the 
friction  coefficient,  C  ,  on  the  bottom  is  assumed  constant  and  if 
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horizontal  mixing  la  neglected,  the  mean  longshore  component  of  velocity, 
,  is  given  as: 

V-  «  (5tt/8)(S/C)  u  sin  a  (9) 

Is  max 

This  value  of  longshore  current  is  proportional  to  the  longshore  com¬ 
ponent  of  the  orbital  velocity.  When  the  horizontal  mixing  is  taken  into 
account  the  longshore  currents  observed  in  prototype  situations  and 
laboratory  experiments  are  consistent  with  a  dimensionless  friction  co¬ 
efficient,  C  ,  of  about  0,01, 

21*  The  profile  of  the  longshore  current,  as  a  function  of 
distance  from  the  swash  line,  was  calculated  by  Longuet-Higgins  (1970a, 
1970b)  by  using  the  concept  of  radiation  stress  together  with  a  hori¬ 
zontal  eddy  viscosity,  ue  ,  of  the  form: 

1/2 

Ue  -  P  N  x  (gd)A/  (10) 

where 

N  -  dimensionless  parameter  that  varied  between  0.000  and  0.016 
x  -  distance  offshore,  ft 

This  development  gives  rise  to  a  family  of  current  profiles  whose  form 
depends  only  on  a  nondimens ional  parameter  H  : 

PL-H  *  77/2  <SN/°*41C)  (ID 

The  current  profiles  are  of  simple  analytic  form,  having  a  maximum  in 
the  surf  zone  and  tending  to  zero  at  the  swash  line.  Comparison  with  the 
laboratory  experiments  of  Galvin  and  Eagleson  (1965)  showed  remarkably 
good  agreement  if  the  friction  coefficient,  C  ,  is  taken  as  0.010.  The 
theoretical  profiles  are  insensitive  to  the  exact  value  of  %  but 

the  experimental  results  suggest  that  P^_^  never  exceeds  a  value  of 
0.016. 

22.  CERC  (1977)  has  modified  the  Longuet-Higgins  (1970a,  1970b) 
expression  for  longshore  velocity,  based  on  more  recent  experimental 


data  and  prototype  observations,  and  determined  that: 


Vj  -  20.7  S  (gHh) 1/2  sin  20^  (12) 

where  is  the  breaker  height.  Here  again  there  is  a  direct  propor¬ 

tionality  to  the  beach  slope,  S  •  Komar  and  Inman  (1970)  concluded 
earlier  that  there  was  no  dependence  on  S  ,  and  proposed  the  relation¬ 
ship  : 

Vls  s  2-7  Um  sin  “b  COS  %  (13) 

in  which 

u  =  (2E. /pd.  ) 1/2  (14) 

m  b  b 

where 

2 

*  wave  energy  density  at  breaking,  ft-lb/ft 
=*  water  depth  at  breaking,  ft 

Komar  and  Inman  (1970)  demonstrated  that  Equation  13  has  fairly  good 
agreement  with  the  available  longshore  current  data,  and  that  it  is 
equivalent  to  the  theoretical  longshore  current  relationship  of  Longuet- 
Higgins  which  was  derived  from  radiation  stress  concepts.  Komar  (1979) 
reiterates  that  the  inclusion  of  the  beach  slope,  S  ,  leads  to  a  some¬ 
what  less  precise  degree  of  fit  for  much  of  the  available  data  (proto¬ 
type  and  laboratory)  than  do  best-fit  expressions  which  do  not  include 
the  beach  slope  parameter. 


Diffraction 

23.  Diffraction  of  water  waves  is  the  phenomenon  by  which  wave 
energy  propagates  into  the  sheltered  lee  of  structures  even  in  the 
absence  of  bathymetric  refraction.  In  these  situations,  wave  crests 
bend  (even  in  constant  depth  water)  and  gradients  of  wave  height 
exist  along  the  wave  crest.  This  phenomenon  is  most  visible  when  a 


train  of  regular  waves  Is  interrupted  by  an  obstruction  such  as  a  jetty 
or  shore-connected  breakwater.  The  theory  of  water  wave  diffraction 
can  be  explained  by  Huygens*  principle.  Each  point  of  an  advancing  wave 
front  (wave  crest)  may  be  considered  as  the  center  of  a  secondary 
circular  wave  which  advances  in  all  directions.  The  resultant  shape  of 
the  crest  is  the  envelope  of  all  these  secondary  waves.  In  a  straight- 
crested  wave  train,  the  envelope  of  the  secondary  waves  is  a  straight 
line  also.  When  the  wave  passes  an  obstruction,  the  energy  intensity 
at  a  certain  point  is  a  vector  combination  of  all  the  circular  waves 
emitted  by  every  point  of  the  passing  wave  train. 

24.  Sommerfeld  (1896)  presented  a  solution  for  the  diffraction  of 
light  waves  past  the  edge  of  a  semi-infinite  screen.  Penny  and  Price 
(1944)  showed  that  this  is  also  the  solution  of  the  water  wave  diffrac¬ 
tion  problem  at  the  end  of  a  semi-infinite  obstacle  such  as  a  jetty  or 
shore-connected  breakwater.  This  exact  solution  of  the  surface  eleva¬ 
tions  behind  the  breakwater  is  applicable  only  to  water  of  constant 
depth  and  waves  of  small  amplitude.  Putman  and  Arthur  (1948)  summarized 
the  solution  of  Penny  and  Price  (1944)  for  the  definitive  sketch  of 
Figure  1.  In  cylindrical  coordinates,  the  water-surface  elevation 
is: 

n  =  (aikc/g)eikct  cosh  kd  F(r,0)  (15) 


where 

r\  *  water-surface  elevation,  ft 
a  *  wave  amplitude,  ft 
c  *  wave  celerity,  ft/sec 
i  =  square  root  of  -1 

F(r,0)  is  a  function  which  satisfies  the  wave  equation  in  cylindrical 
coordinates: 
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Figure  1.  Definitive  sketch,  wave  diffraction 
around  a  breakwater 


In  the  presence  of  a  jetty  or  breakwater,  the  boundary  condition  is 
imposed  that  the  normal  component  of  the  fluid  velocity  is  zero  along 
the  breakwater,  leading  to  the  solution: 
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Bretschneider  (1966)  has  presented  computational  procedures  for  evalu¬ 
ating  the  diffraction  coefficients  at  arbitrary  points  behind  jetties 
or  breakwaters • 

25.  Wiegel  (1962)  developed  a  graphical  procedure  for  determining 
diffraction  coefficients  of  waves  passing  the  tip  of  single  breakwaters. 
The  family  of  diagrams  shows,  for  uniform  water  depth,  lines  of  equal 
wave-height  reduction  displayed  in  terras  of  the  diffraction  coefficients. 
The  diffraction  diagrams  (typical  example,  Figure  2)  are  constructed  in 
polar  coordinate  form  centered  at  the  structure  tip.  The  arcs  behind 
the  breakwater  are  spaced  one  radius-wavelength  unit  apart  so  that, 
in  application,  a  specific  diagram  must  be  scaled  up  or  down  so  that 
the  particular  wavelength  corresponds  to  the  scale  of  the  hydrographic 
area  under  investigation.  The  set  of  diffraction  diagrams  of  waves 
passing  the  tip  of  a  single  breakwater  were  presented  by  CERC  (1977), 
along  with  diagrams  of  waves  entering  a  breakwater  (typical  example, 
Figure  3). 

26*  Laboratory  tests  were  performed  by  Harms  (1979)  to  investigate 
the  distribution  of  wave  heights  in  the  lee  of  a  breakwater  (shore- 
connected)  for  waves  normally  incident  upon  the  structure  and  with  a 
horizontal  bottom  both  in  front  of  and  in  the  lee  of  the  structure.  In 
general,  satisfactory  agreement  was  obtained  between  measurement  and 
theory,  but  diffraction  theory  was  not  found  to  be  conservative.  At 
large  distances  in  the  shadow  zone,  measured  wave  heights  consistently 
exceeded  theoretical  values.  Close  to  the  breakwater  outside  the  shadow 
zone,  the  measured  maximum  wave  height  was  also  found  to  be  larger  than 
that  predicted  by  theory.  The  diffraction  behavior  appeared  to  be  in¬ 
sensitive  to  the  intensity  of  wave  reflections  from  the  seaward  side  of 
the  breakwater. 
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Figure  3.  Generalized  diffraction  diagram  for  a  brea 

wavelengths  (after  CERC  1977) 


Combined  Refraction  and  Diffraction 


27.  The  bathymetry  shoreward  of  a  jetty  or  breakwater  usually  is 
not  flat  or  even  uniform;  hence,  refraction  generally  occurs  in  addition 
to  the  diffraction  effects.  While  a  general  unified  analytical  approach 
to  the  simultaneous  solution  of  these  two  distinctly  different  phenomena 
has  not  been  entirely  developed,  considerable  insight  has  been  gained 
through  the  theoretical  work  of  Liu  and  Mei  (1975,  1976a,  1976b),  Lick 
(1978),  Liu  and  Lozano  (1979),  and  through  the  earlier  experimental  work 
of  Mobarek  (1962).  The  procedure  usually  followed  by  coastal  engineers 
concerned  with  wave-height  variation  behind  jetties  or  breakwaters  is  to 
construct  refraction  diagrams  shoreward  to  the  structure,  then  construct 
diffraction  diagrams  for  three  or  four  wavelengths  shoreward  of  the 
jetty,  and  finally  refract  the  last  wave  crest  on  toward  the  shoreline. 
This  procedure  is  schematized  in  Figure  4,  where  the  overall  refraction 
diffraction  coefficient,  ^  ,  in  the  region  behind  the  structure  is: 


where 

K  *  refraction  coefficient  at  the  structure,  dimensionless 
r 

*  diffraction  coefficient  on  last  wave  crest  behind  the 

structure  from  which  additional  refraction  computations,  are 
performed,  dimensionless 

bj  »  orthogonal  spacing  at  the  last  diffracted  wave  crest,  ft 
b2  *  orthogonal  spacing  near  the  shore,  ft 

Wave  heights  resulting  from 
refraction  and  diffraction 

28.  Mobarek  (1962)  experimentally  investigated  the  effect  of 
bottom  slope  on  wave  diffraction  through  a  gap  in  a  breakwater  normal  to 
the  incident  wave  direction.  Also  Investigated  was  the  effect  of  an 
abrupt  increase  or  decrease  in  the  water  depth  behind  the  breakwater. 

The  theoretical  analysis  for  the  comparison  of  experimental  results  fol¬ 
lowed  the  treatment  of  Penny  and  Price  (1944)  restricted  to  the  case  of 


Figure  4,  Graphical  procedure  for  determining  overall  refraction  and 
diffraction  effects  past  a  breakwater  (after  CERC  1977) 

normal  incidence  for  which  the  Sommerfeld  (1896)  solution  is  simplified 
and  in  the  presence  of  a  horizontal  bottom.  Two  fundamentally  different 
basin  configurations  were  used  in  the  study.  The  first  (Figure  5)  con¬ 
sists  of  a  longitudinally  sloping  bottom  with  the  slope  commencing  at 
the  breakwater  and  extending  to  the  shoreline.  The  second  (Figure  6) 
was  constructed  with  a  flat  bottom  extending  beyond  the  breakwater  gap 
but  sloping  laterally  to  the  shoreline. 

29.  Results  of  the  Mobarek  (1962)  investigation  for  a  longitudi¬ 
nally  sloping  bottom  and  for  a  laterally  sloping  bottom  are  shown  in 
Figures  7  and  8,  respectively.  Taking  into  consideration  the  serious 
limitations  of  the  experimental  equipment  (very  small  model,  72  sq  ft), 
the  investigation  led  to  the  conclusion  that  the  procedure  usually 
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Figure  7.  Combined  refraction  and  diffraction  coefficients,  Kr_<j  , 
for  longitudinal  bed  slope  behind  breakwater  (after  Mobarek  1962) 
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followed  for  estimating  wave  heights  behind  jetties  or  shore-connected 
breakwaters  (Figure  4)  was  sufficiently  good  for  medium  period  waves; 
but  in  the  case  of  long  waves,  the  effect  of  the  shoaling  bottom  on 
waves  should  be  taken  into  consideration*  Experiments  on  a  larger 
scale  were  highly  recommended* 

30.  For  a  long  shore-connected  breakwater  on  a  slowly  varying 
bottom,  an  asymptotic  theory  has  been  developed  by  Liu  and  Mei  (1976a) 
which  accounts  for  the  combined  effects  of  refraction  and  Fresnel  dif¬ 
fraction  of  water  waves.  Qualitatively,  the  axis  of  the  diffraction 
zone  downwave  of  the  breakwater  or  jetty  was  found  to  be  curved;  and 
quantitatively ,  the  diffracted  wave  amplitude  must  be  modulated  hori¬ 
zontally  an  amount  depending  on  the  refracted  wave  amplitude*  The  total 
wavefield  is  complex,  especially  near  the  tips  of  the  jetty,  and  is 
caused  by  the  interference  between  the  incident  and  the  reflected  waves. 
An  example  solution  is  presented  in  Figure  9  for  a  jetty  assumed  to  be 
400  m  long  from  the  still-water  shoreline.  The  beach  slope  was  0.1  for 
x  <  100  m  and  0.01  for  x  >  100  m  •  The  incident  wave  period  was 
10  sec,  and  the  initial  amplitude  was  0.5  m  with  the  incidence  at  in¬ 
finity  of  45  deg.  For  short  jetties  or  groins,  the  reflection  and  dif¬ 
fraction  effects  may  be  either  too  localized  or  too  complicated  by 

METRES 


Figure  9.  Variation  of  the  diffraction  factor  downwave  of  a 
shore-connected  breakwater,  incidence  angle  ■  45  deg  (after 

Liu  and  Mei  1976a) 
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virtue  of  the  dominant  influence  of  the  surf  zone  to  be  handled  by  this 
theory.  The  approximation  is  a  combination  of  Kirchhoff's  theory  and 
geometrical  optics.  This  permits  the  diffraction  field  near  the  beach 
to  be  calculated  analytically  and  this  can,  in  turn,  be  used  for  breaker 
line  estimates. 

31.  A  uniformly  valid  asymptotic  solution  for  water  waves  has 
been  developed  by  Liu  and  Lozano  (1979)  that  accounts  for  the  combined 
effects  of  refraction  due  to  slowly  varying  water  depth  and  diffraction 
by  a  long  shore-connected  breakwater.  This  solution  is  more  general 
than  the  approximate  solution  developed  by  Liu  and  Mei  (1976a)  because 
this  theory  is  valid  near  the  edge  and  the  tip  of  the  breakwater.  The 
wave  behavior  in  the  near  field  is  of  particular  interest  for  studying 
the  scour  and  erosion  that  may  occur  near  the  tip  of  a  breakwater.  In 
this  analysis,  recent  developments  in  the  field  of  formal  geometrical 
optics  have  been  extended  to  include  the  effects  of  diffraction.  A 
numerical  example  is  presented  for  a  linear  sloping  bottom,  although 
the  theory  is  not  restricted  to  this  condition.  The  length  of  the  break¬ 
water  is  assumed  to  be  400  m  positioned  on  a  beach  slope  of  IV  on  20H. 

The  incident  wave  train  with  period  of  10  sec  arrives  at  an  angle  of 
45  deg.  Figure  10  shows  dimensionless  results  of  the  wave  amplitude 
variation  downwave  of  the  breakwater  for  two  sections  taken  perpendicular 
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Figure  10.  Dimensionless 
wave  amplitudes  perpendic- 
cular  to,  and  downwave  of, 
shore-connected  breakwater 
(after  Liu  and  Lozano  1979) 
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to  the  breakwater  at  distances  of  100  and  200  m  from  the  shoreline. 

32.  Berkoff  (1972)*  Schonfeld  (1972),  and  Smith  and  Sprinks  (1975) 
independently  derived  a  two-dimensional  equation  that  governs  short-wave 
propagation  over  moderately  varying  depths.  Rewriting  Berkoff's  expres¬ 
sion,  Jonsson  and  Brink-Kjaer  (1973)  Introduced  the  mild-slope  wave 
equation: 
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*  horizontal  gradient  operator,  dimensionless 
=  phase  velocity,  ft/sec 

=  group  velocity,  ft/sec 

*  angular  frequency,  2tt/L,  1/ft 
=  complex  wave  amplitude,  ft 

21  can  be  rewritten  in  the  form: 
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where  <t>  is  a  velocity  potential  defined  by  u  »  V<}>  and  u  is  a  two- 
dimensional  vector.  Equation  22  reduces  to  the  diffraction  Helmholtz 
equation  in  deep  or  constant-depth  water.  In  shallow  water,  the  equation 
reduces  to  the  linear  long-wave  equation. 

33.  Houston  (1980)  solved  Equation  22  by  the  use  of  a  hybrid 
finite-element  numerical  model  originally  developed  by  Chen  and  Mei 
(1974)  to  solve  the  diffraction  Helmholtz  equation  in  a  constant-depth 
region.  The  appropriate  modifications,  including  variable  depth  and 
frequency  dispersion,  were  Incorporated  by  Houston  (1980)  and  the  solu¬ 
tion  of  Equation  22  was  applied  to  the  geometry  of  the  physical  model 
(described  in  PARTS  III  and  IV  of  this  report)  used  in  this  investiga¬ 
tion.  A  problem  in  simulating  these  hydraulic  model  tests  numerically 
is  that  the  waves  break  in  the  hydraulic  model  near  the  shoreline  and 
thus  dissipate  their  energy.  No  mechanism  exists  to  dissipate  energy 
in  the  numerical  model.  However,  dissipation  was  simulated  by  allowing 
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waves  to  continue  to  propagate  out  of  the  problem  area.  The  breakwater 
and  uniform  slope  were  modeled  only  to  the  point  where  breaking  occurred. 
The  depth  was  then  increased  to  the  depth  of  a  semi-infinite  ocean 
region  surrounding  the  region  of  computation,  and  the  waves  were  allowed 
to  radiate  away  from  the  area  of  interest. 

34.  Figure  11  shows  a  typical  comparison  between  the  hydraulic 


DISTANCE  NORMAL  TO  BREAKWATER,  FT 

Figure  11.  Comparison  of  asymptotic  and  finite- 
element  solutions  with  hydraulic  model  measurements 
(after  Houston  1980) 

model  measurement's  and  the  finite-element  calculations  (Houston  1980). 
Also  shown  is  the  uniformly  valid  asymptotic  solution  derived  by  Liu  and 
Lozano  (1979).  The  solution  derived  by  Liu  and  Lozano  (1979)  appears  to 
be  in  good  agreement  with  the  hydraulic  model  tests.  The  finite-element 
calculations  of  Houston  (1980)  agree  quite  well  in  the  shadow  zone  with 
the  hydraulic  model  tests.  The  agreement  is  not  as  good  outside  the 
shadow  zone.  The  difference  is  partially  attributable  to  the  artificial 
increase  in  depth  to  allow  the  waves  to  radiate  from  the  inner  region. 
This  depth  transition  would  cause  some  energy  to  reflect  back  into  the 
region  of  interest. 


Wave-induced  currents  produced 
by  refraction  and  diffraction 


35.  The  wave  field  near  a  shore-connected  breakwater  has  been 
studied  by  Liu  and  Mei  (1975,  1976b)  for  the  purpose  of  developing  the 
wave-induced  current  patterns,  magnitudes,  and  circulation  cells  in  the 
neighborhood  of  the  breakwater.  By  employing  an  approximation  which 
adequately  describes  Fresnel  diffraction,  the  first  order  wave  was  de¬ 
duced  accounting  for  diffraction  and  refraction.  In  a  natural  environ¬ 
ment,  the  waves  and  current  change  the  bottom  topography  and  vice  versa, 
with  the  actual  beach  process  being  one  of  interaction.  Efforts  were 
focused  on  the  effect  of  wave  pattern  due  to  reflection,  refraction,  and 
diffraction  on  the  current,  assuming  the  topography  to  be  fixed  and 
simple.  The  beach  was  considered  to  be  immobile,  infinitely  long,  and 
uniform  in  the  longshore  direction.  The  depth  contours  were  straight 
and  parallel.  The  walls  of  the  breakwater  were  vertical,  impermeable, 
perfectly  reflecting,  and  of  negligible  thickness  compared  with  a  wave¬ 
length  and  extended  far  into  the  surf  zone.  For  mathematical  simplifica¬ 
tion,  the  assumptions  were  made  that  turbulent  lateral  stresses  and 
convective  inertia  terras  could  be  ignored.  The  bottom  friction  was  as¬ 
sumed  proportional  to  the  square  of  the  average  current  velocity. 

36.  The  assumption  by  Liu  and  Mei  (1975)  of  the  law  of  bottom 
friction  has  been  made  by  other  researchers  and  the  constant  of  propor¬ 
tionality  (friction  factor)  determined  by  laboratory  and  field  experi¬ 
ments.  Longuet-Higgins  (1970a)  and  Thornton  (1970)  investigated  typical 
ranges  of  possible  values  of  the  friction  factor  based  on  laboratory  ex¬ 
periments  with  rigid  and  rough  bottoms.  Typical  values  appeared  to  lie 
within  the  range  0.005  to  0.012.  In  field  conditions,  other  factors 
such  as  sediment  motion  at  the  bottom  may  enter  into  consideration. 
Bretschneider  (1954)  reported  from  measurements  in  shallow  water  in  the 
Gulf  of  Mexico  that  for  wind  waves  the  friction  factor  varied  from  0.03 
to  0.097,  with  an  average  value  of  0.053.  Based  on  laboratory  and  field 
studies,  Putnam,  Munk,  and  Traylor  (1949)  reported  values  between  0.07 
and  0.385.  For  comparison  purposes,  Liu  and  Mei  (1975)  used  two  values, 
0.01  and  0.09. 
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37.  Sample  calculations  were  prepared  by  Liu  and  Mei  (1975),  but 
because  of  the  many  scales  involved,  the  results  are  presented  in  physi¬ 
cal  dimensions.  A  shore-connected  breakwater  of  400  m  length  was  evalu¬ 
ated,  with  the  incident  wave  being  of  10-sec  period  and  having  an  ampli¬ 
tude  of  0.5  m.  The  angle  of  incidence  in  deep  water  was  45  deg.  One  of 
the  more  significant  results  of  the  study  was  the  formation  of  cells. 
Whether  in  the  diffraction  zone  (downwave)  or  in  the  reflection  zone, 
these  cells  are  associated  with  variations  in  the  breaking  wave  height 
which  is  closely  related  to  the  variations  in  the  breaker  line  position, 
the  surf  zone  width,  and  the  mean  sea  level.  In  the  downwave  zone,  the 
current  generally  follows  the  breakwater  and  the  shore.  Around  the  point 
where  the  breaking  wave  is  the  highest,  there  is  a  large  counterclockwise 
cell  near  the  shore.  The  corresponding  mean  sea  level  shows  that  little 
disturbance  is  felt  near  the  breakwater.  Near  the  breaker  line,  the  mean 
sea  level  is  lowest.  As  the  longshore  current  is  approached,  the  mean  sea 
level  contours  become  parallel  to  the  shoreline.  The  streamline  patterns 
for  coefficients  of  friction  of  0.01  and  0.09  are  shown  in  Figures  12 


Figure  12.  Streamline  pattern,  downwave  region,  C  -  0.01 
(after  Liu  and  Mei  1975) 
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Figure  13.  Streamline  pattern,  downwave  region,  C  «  0.09 
(after  Liu  and  Mei  1975) 

and  13,  respectively,  for  the  downwave  region.  Contours  of  mean  sea 
level  for  the  downwave  region  are  shown  in  Figure  14. 

38.  Liu  and  Mei  (1975)  were  not  able  to  find  any  direct  observa¬ 
tion  of  currents  comparable  to  the  shore-connected  breakwater  problem 
which  they  had  investigated  analytically.  The  existing  experiments  at 
that  time  dealt  largely  with  configurations  of  a  harbor  breakwater  or  a 
pair  of  inlet  jetties.  Short  groins  were  not  strictly  relevant  to  the 
study  of  Liu  and  Mei  (1975).  For  short  groins,  the  diffraction  effect 
could  be  overwhelmed  by  the  inertial  effect  of  the  incoming  longshore 
current  which  would  be  simply  deflected  seaward  as  a  jet.  Hence,  cells 
near  the  shore  should  not  be  expected  to  form.  To  provide  insight  and 
gain  a  better  understanding  of  the  phenomena,  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  performed  physical  model  investiga¬ 
tions  of  a  shore-connected  single  breakwater  oriented  perpendicular  to 
the  shoreline.  Wave-height  distributions,  current  magnitudes  and  pat¬ 
terns,  and  cell  formations  were  determined  and  are  described  in  PART  IV. 
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PART  III:  EXPERIMENTAL  FACILITIES 


Physical  Model 

39-  This  experimental  investigation  was  conducted  in  the  fa¬ 
cilities  of  the  Hydraulics  Laboratory  of  WES.  The  physical  model, 
which  was  molded  in  cement  mortar,  consisted  of  a  50-  60-ft  area, 

with  a  water  depth  of  1  ft  in  the  open-ocean  region  (Figure  15).  The 
beach  slope  was  IV  on  20H,  with  a  vertical,  impermeable  breakwater 
extending  perpendicularly  from  the  shoreline  for  a  distance  of  15  ft 
to  a  point  where  the  water  depth  was  0-75  ft.  The  uniform  slope 
continued  on  to  a  water  depth  of  1  ft,  beyond  which  the  basin  was 
horizontal  to  the  facility  walls.  The  wave  generator  was  mobile, 
and  two  positions  were  tested  (incident  waves  of  20  and  30  deg) • 

Wave  guides  extended  from  the  ends  of  the  generator  along  wave  or* 
thogonals  toward  the  shoreline-  Orthogonals  were  determined  by 
standard  refraction  techniques  for  waves  propagating  onto  a  uniform 
slope  in  the  absence  of  the  breakwater.  The  area  of  interest  was 
downwave  of  the  breakwater-  Four  sections  perpendicular  to  the  break¬ 
water  were  instrumented  with  wave-height  sensors  for  determination 
of  the  wave  field-  These  sections  were  located  6,  8,  10,  and  12  ft 
from  the  mean  waterline  (shoreline  end  of  the  breakwater),  and  are 
labeled  "distance  along  the  breakwater-" 

40.  Model  waves  were  generated  by  a  30-ft-long  wave  generator 
with  a  trapezoidal-shaped,  vertical-motion  plunger-  The  vertical 
movement  of  the  plunger  caused  a  periodic  displacement  of  water 
surface.  The  length  of  stroke  and  the  frequency  of  the  vertical 
motion  were  variable  over  the  range  necessary  to  generate  waves 
with  the  required  characteristics.  In  addition,  the  wave  generator 
was  mounted  on  retractable  casters  which  enabled  it  to  be  positioned 
to  generate  waves  from  the  required  directions.  A  2-ft  (horizontal) 
solid  layer  of  fiber  wave  absorber  was  placed  around  the  inside 
perimeter  of  the  model  to  damp  wave  energy  that  might  otherwise 
be  reflected  from  the  model  walls. 
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Figure  15.  Experimental  facility 


Instrumentation  for  Wave-Height  Determination 

41.  The  increase  in  recent  years  in  the  number  and  complexity  of 
wave  model  studies  conducted  by  the  Hydraulics  Laboratory  of  WES  and 
the  use  of  such  models  to  solve  long-wave  problems  (requiring  large  model 
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areas)  have  demonstrated  a  need  for  automation  of  modeling  procedures 
including  operation,  data  collection,  and  data  analysis  (Durham  and 
Greer  1976).  Automated  Data  Acquisition  and  Control  Systems  (ADACS) 
have  been  designed  and  built  at  WES  for  the  purpose  of  collecting  wave 
data  and  controlling  operations  of  hydraulic  wave  models.  The  computer 
hardware  configuration  for  each  system  consists  of  a  minicomputer  with 
32k  16-bit  words  of  memory,  a  magnetic  tape  controller  with  two  9-track 
tape  drives,  one  moving  head  disc  controller  with  one  removable  platter 
and  one  nonremovable  platter,  an  interval  timer  (1  nsec) ,  an  analog  to 
digital  12-bit  converter  featuring  64  analog  (+10  volts)  inputs  and  a 
45  kHz  multiplexer,  an  ASR  33  teletype  unit,  96  sense/control  lines,  and 
one  matrix  electrostatic  printer/plotter. 

42.  The  ADACS  are  capable  of  automatically  calibrating  the  wave 

sensors,  controlling  wave  generators,  acquiring  data  from  the  sensors  at 
a  high  sampling  rate,  and  analyzing  test  data.  Data  are  taken  and  re¬ 
corded  on  disc  or  magnetic  tape  for  direct  analyses  by  the  minicomputer 
system  or  on  magnetic  tape  in  a  format  compatible  with  a  Honeywell  G635 
for  backup  analyses.  Automatic  calibration  of  wave  sensors  has  reduced 
the  time  required  to  calibrate  the  sensors  by  a  factor  of  four.  In 
addition,  several  times  the  number  of  tests  can  be  run  during  a  day  with 
test  results  analyzed  at  completion  of  model  tests  by  minicomputer.  The 
system  configuration  (Figure  16)  of  ADACS  consists  of  the  following  sub¬ 
systems:  (a)  digital  data  recording  and  controls,  (b)  analog  recorders 

and  channel  selection  circuits,  (c)  wave  sensors  and  interfacing  equip¬ 
ment,  and  (d)  wave  generators  and  control  equipment. 

Wave  sensors 

43.  The  data  acquired  from  wave  models  are  the  water-surface 
variations  about  a  reference  water  level.  This  information  is  col¬ 
lected  at  selected  geographic  locations  within  the  model  for  specified 
wave  conditions  at  the  wave  generator.  Wave  sensors  are  used  to  obtain 
this  information  at  selected  locations  in  the  model.  Each  of  the 
water-surface-piercing,  parallel  rod  wave  sensors  is  connected  to  a 
Wheatstone  bridge  and  a  transducer  measures  the  conductance  of  the 
water  between  the  two  parallel  rods  which  are  mounted  vertically 
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Figure  16.  Schematic  of  components  of  Automated  Data 
Acquisition  and  Control  System  (ADACS) 


(Figure  17).  The  conductance  is  directly  proportional  to  the  depth  of 
submergence  of  the  two  rods  in  the  water.  The  output  of  each  wave  sensor, 
is  routed  through  shielded  cables  to  its  signal  conditioning  equipment 
where  it  is  processed  for  recording.  The  ADACS  can  detect  changes  in 
water-surface  elevations  to  an  accuracy  of  0.001  ft.  To  obtain  this 
accuracy*  ultrastable  power  supplies  and  good  signal-to-noise  ratios  are 
necessary.  The  carrier  source  for  the  wave  sensor  bridge  maintains  a 
variation  of  less  than  0.025  percent. 

Calibration 


44. 


In  order  to  convert  the  water-elevation  data  in  millivolts  to 
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Figure  17,  Parallel  rod  conductance-type  wave  sensor 

water-surface  elevations  in  feet,  each  wave  sensor  must  be  calibrated. 

The  capability  of  automatically  calibrating  the  wave  sensors  (maximum  of 
25  rods  simultaneously)  prior  to  collecting  data  is  provided  by  ADACS. 

To  calibrate  each  set  of  parallel  rods,  the  voltage  from  the  signal  con¬ 
ditioning  equipment  is  monitored  and  recorded  as  the  parallel  rods  are 
moved  vertically  a  known  distance  into  or  out  of  the  water.  A  precision, 
linear-position  potentiometer  is  located  on  the  wave  sensor  stand  and  is 
coupled  directly  to  the  parallel  rods  by  a  gear  train  driven  by  electric 


motor.  By  moving  vertically  the  coupled  wave  sensor  and  potentiometer 
with  the  electric  motor  and  by  monitoring  the  output  voltage  from  the 
potentiometer,  the  wave  sensor  can  be  moved  vertically  a  precise 
distance.  The  electric  motor  for  each  wave  sensor  is  controlled  by 
a  control/sense  line  and  a  relay  contact.  The  minicomputer  controls  the 
vertical  movement  of  each  wave  sensor  by  actuating  the  control /sense 
line.  The  central  processing  unit  acts  as  a  voltage  comparator  by  moni¬ 
toring  the  potentiometer  voltage  and  comparing  it  with  a  reference 
voltage  which  is  determined  from  desired  displacement  and  potentiometer 
calibration.  When  the  voltage  comparison  is  satisfied,  the  control/ 
sense  line  is  reactivated,  the  electric  motor  stops,  and  voltage 
samples  from  the  rods  and  potentiometers  are  acquired.  By  systemati¬ 
cally  moving  each  wave  sensor  through  11  quasi-equal ly  spaced  locations 
over  the  range  of  rod  length  used,  voltage  versus  known  displacements 
are  obtained  from  which  a  calibration  curve  for  each  sensor  can  be 
calculated  and  recorded  on  magnetic  tape  or  disc.  After  collecting 
the  calibration  data,  the  minicomputer  analyzes  these  data  by  least- 
squares  fitting  a  set  of  curves  (linear,  quadratic,  or  spline)  to  the 
data,  determining  the  best  order  of  fit,  and  comparing  the  maximum 
deviation  of  the  best  fit  with  a  previously  acceptable  value  for  this 
maximum  deviation. 

Data  acquisition  and  analysis 

45.  During  the  acquisition  mode,  wave  data  for  a  specified  wave 
condition  at  the  wave  generator  are  collected  from  a  maximum  of  50  wave 
sensors,  recorded  on  analog  strip  charts,  digitized,  and  recorded  on 
magnetic  tape  or  disc  for  further  analyses.  The  sampling  scheme  is 
flexible  and  can  be  tailored  for  different  applications  with  maximum 
throughput  rates  theoretically  limited  by  the  multiplexer  rate  and 
allocatable  buffer  size.  The  sampling  scheme  used  in  this  investigation 
was  60  discrete  voltage  samples  equally  spaced  over  each  wave  period  for 
a  predetermined  number  of  90  wave  periods,  for  each  of  the  sensor  loca¬ 
tions.  The  minicomputer  calculated  from  Input  parameters  the  lag  at  the 
beginning  of  data  acquisition  by  10  wave  periods  after  starting  the 
generator,  provided  timing  pulses  for  synchronizing  and  controlling  the 
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recorders*  and  determined  completion  of  the  test*  The  determination  of 
the  height  of  each  wave  of  the  monochromatic  wave  train  was  performed 
(at  each  sensor  location),  the  average  of  these  90  individual  heights 
was  calculated*  and  the  standard  deviation  of  these  individual  observa¬ 
tions  about  the  mean  was  computed.  The  value  displayed  as  the  wave 
height  at  each  sensor  location  is  this  mean  value  plus  or  minus  one 
standard  deviation.  The  data  were  obtained  by  operating  32  wave-height 
sensors  simultaneously  for  each  test  condition*  and  then  repeating  the 
same  test  with  the  sensors  repositioned  to  allow  better  definition  of 
wave-height  gradients  away  from  the  breakwater.  A  typical  arrangement 
of  the  32  sensors  near  the  breakwater  is  shown  in  Figure  18 •  Only  the 
wave  sensor  probe  penetrates  the  water  surface*  thereby  minimizing  local 
disturbances. 


Figure  18.  Wave-height  sensor  arrangement,  grid  spacing  *  2  ft 
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Average  Current  Determination 


46,  To  define  the  circulation  patterns  downwave  of  the  breakwater v 
average  values  of  wave-induced  velocity  were  determined  at  those  grid 
locations  shown  in  Figure  19*  The  average  values  were  determined  at 
each  individual  station  after  the  wave  generator  had  been  operating  for 

2  min  so  that  the  circulation  cell  immediately  adjacent  to  the  break¬ 
water  could  become  well  established  prior  to  measurements.  After  the 
measurement  at  either  surface,  middepth,  or  bottom  had  been  obtained  at 
each  station,  wave  generator  operation  was  terminated  so  that  those 
fluid  motions  unique  to  the  physical  model  (downcoast  circulation  cell 
Induced  by  the  boundary  of  the  finite  model  area)  would  not  have  time  to 
progress  into  the  area  of  major  interest  (the  immediate  vicinity  of  the 
structure).  It  was  known  a  priori  that  the  existence  of  the  downcoast 
boundary  condition  would  preclude  exact  duplication  of  an  infinite  coast¬ 
line;  however,  for  the  purpose  of  verifying  various  theories,  numerical 
models  can  be  formulated  to  represent  a  finite  model  area  that  Includes 
all  boundaries  appropriately  defined. 

47.  Surface  current  magnitudes  were  determined  by  timing  a 
particle  floating  on  the  surface  and  initially  positioned  over  the  grid 
points.  The  time  and  direction  required  for  each  particle  to  traverse 
a  distance  of  1  ft  was  noted  and  recorded.  Surface  velocities  were 
especially  sensitive  to  wave  orbital  motions.  Nlddepth  velocities  were 
determined  by  observing  the  movement  of  dye  placed  at  one-half  the  still- 
water  depth.  While  a  measure  of  dispersion  would  tend  to  occur  after  a 
finite  time  interval,  the  center  of  the  dye  region  could  be  easily  ascer¬ 
tained;  and  this  was  the  element  considered  in  the  average  velocity  de¬ 
terminations.  Bottom  current  velocities  and  directions  were  determined 
by  observing  dye  movement  which  had  been  placed  essentially  on  the  bottom 
of  the  model.  Except  for  those  regions  in  the  breaker  rone,  the  dye 
patterns  on  the  bottom  were  more  easily  observable  than  those  at  mid- 
depth  because  downward  dispersions  were  precluded.  Hence,  the  dye  region 
moved  in  a  more  concentrated  form. 
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PART  IV:  RESULTS 


Wave  Heights 

48.  The  experimental  facility  used  in  this  study  was  constructed 
so  that  transition-to-shal low-water  waves  would  exist  near  the  breakwater 
(representative  of  prototype  waves  which  are  capable  of  causing  scour 
and  erosion).  This  requirement  indicated  that  transition-depth  waves 
should  be  generated,  since  physical  model  size  and  depth  constraints  pre¬ 
cluded  the  generation  of  deepwater  waves  with  required  characteristics. 

In  addition,  the  requirement  existed  that  the  generated  wave  heights 
should  be  large  enough  so  that  small  changes  in  the  heights  could  be  de¬ 
tectable,  and  at  the  same  time  the  heights  should  be  small  enough  so 
that  comparisons  with  linear  theories  could  be  performed.  Preliminary 
tests  indicated  that  for  the  area  of  major  interest  and  for  the  range  of 
wave  periods  considered  pertinent,  the  specific  test  conditions  shown  in 
Table  1  could  be  experimentally  investigated  with  height  changes  remain¬ 
ing  essentially  linear,  thus  permitting  comparisons  with  theoretical 
developments. 

49.  Typical  examples  of  the  manner  in  which  the  wave  profile 
transforms  under  the  combined  effects  of  refraction  and  diffraction  are 
displayed  In  Figure  20  for  a  1.00-sec  wave  approaching  the  breakwater 
from  an  incident  direction  of  20  deg.  Positions  of  the  two  referenced 
gages  are  shown  in  Figure  21.  Typical  wave  patterns  covering  the  range 
of  conditions  tested  from  the  20  deg  incident  direction  are  shown  in 
Photos  1-15.  Wave-height  amplification  coefficients,  H/Hq  *  for  all 
conditions  tested,  are  shown  in  Plates  1-48.  Each  plate  constitutes  a 
section  perpendicular  to  the  breakwater  at  distances  from  the  shoreline 
of  6,  8,  10,  or  12  ft,  and  depicts  the  average  maximum  wave  height  that 
occurred  at  each  sensor  location.  The  corresponding  scatter  of  the 
experimental  data  is  Indicated  by  the  error  bars  of  plus-or-minus  one 
standard  deviation. 

50.  Thirty-four  wave-height  sensors  were  used  to  determine  the 
wave  heights  along  the  four  sections  perpendicular  to,  and  downwave  of. 


,7  ,*-*</  •yrAvW’i 


46 


Figure  20.  Typical  wave  form  transformation  after  leaving  generator,  wave  period  *  1.00  sec 
wave  height  near  generator  *  0.112  ft,  approach  angle  of  incidence  -  20  deg 


Figure  21.  Location  of  gages  6  and  14 


the  vertical  breakwater  (eight  gages  along  each  section  and  two  reference 
gages  in  the  ocean  region  near  the  wave  generator) .  The  average  of  the 
heights  recorded  at  the  two  ocean  gages  was  selected  to  be  the  input 
wave  height  from  the  wave  generator.  The  heights  recorded  at  the  re¬ 
maining  32  gages  were  normalized  to  this  input  wave  height.  All  34  gages 
were  recorded  simultaneously  with  each  individual  reading  at  all  data 
stations  consisting  of  the  average  of  90  waves.  Repeatability  testing 
was  conducted  under  identical  test  conditions  to  define  the  variability 
of  the  measurements.  The  numbers  of  replications  varied  from  3  to  10 
and  are  shown  in  Table  2.  The  statistical  measure  of  the  variability 
was  the  square  root  of  the  variance,  or  the  standard  deviation.  The 
variance  is  defined  as  the  sum  of  the  squares  of  the  deviates  of  each 
individual  observation  from  their  average,  divided  by  one  less  than  the 
total  number  of  deviates.  The  average  of  the  ratios  is  displayed  in 
Plates  1-48,  with  error  bars  showing  plus  or  minus  one  standard  devia¬ 
tion.  One  standard  deviation  is  usually  2  to  3  percent  of  the  average 
value  of  the  observations.  The  data  were  originally  obtained  with  the 
gages  positioned  at  2-ft  intervals  starting  at  the  breakwater  and  ex¬ 
tending  to  a  location  14  ft  from  the  breakwater  (even  number  of  feet 
from  the  breakwater  in  Table  2).  To  provide  better  definition  of  the 
wave  field,  the  gages  were  repositioned  at  2-ft  intervals  starting  at  a 
location  1  ft  from  the  breakwater  and  extending  to  a  location  15  ft  from 
the  structure  (odd  number  of  feet  from  the  breakwater  in  Table  2).  Also 
shown  in  Plates  1-20  are  the  wave  heights  determined  by  diffraction 
theory  alone,  based  on  the  semigraphica 1  constant  depth  procedure  of  CERC 
(1977).  These  were  the  only  identical  situations  available  for  compari¬ 
son.  The  fit  is  seen  to  be  best  at  most  locations  for  short-period 
waves.  Agreement  also  is  best  at  the  seaward  end  of  the  breakwater, 
since  refraction  becomes  more  significant  the  farther  the  waves  propa¬ 
gate  toward  the  shoreline. 

51.  The  incident  wave  heights  generated  in  the  model  were  suf¬ 
ficiently  small  such  that  the  waves  were  usually  linear  within  the  mea¬ 
surement  region  shown  in  Figure  15.  The  degree  of  linearity  is  displayed 
in  Tables  3-14  for  representative  replications  of  each  different  wave 
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characteristic  (period,  height,  and  direction  of  approach)  tested.  It  was 
found  that  at  the  ocean  gages  on  the  average,  98.50  percent  of  the  energy 
of  the  wave  form  exists  at  the  generated  period  (fundamental  frequency). 

Of  the  remaining  1.50  percent  of  the  energy,  1.28  percent  exists  at  the 
first  harmonic.  The  distribution  of  the  wave  energy  throughout  the  model 
for  these  representative  experiments  is  displayed  in  Table  3-14  also. 

Here  is  seen  the  manner  in  which  nonlinear  effects  are  detectable  as  the 
wave  propagates  shoreward  past  the  various  gages.  A  portion  of  the 
energy  of  the  wave  form  is  redistributed  from  the  fundamental  frequency 
to  higher  harmonics,  and  only  95.25  percent  of  the  energy  exists  (on  the 
average)  at  the  generated  period.  However,  of  the  remaining  4.75  percent 
of  the  energy,  4.16  percent  exists  (on  the  average)  at  the  first  harmonic. 

52.  In  the  neighborhood  of  the  breakwater,  there  were  currents 
(described  in  the  next  section)  that  affected  the  wave  heights.  Oppos¬ 
ing  currents  are  known  to  produce  significant  increases  in  wave  heights 
for  relatively  small  values  of  current-to-wave  velocity  ratios.  Since 
the  magnitude  <  f  the  currents  is  a  function  of  the  incident  wave  height, 
the  ratio  of  measured  wave  height  to  incident  wave  height,  H/H^  >  is 
expected  to  be  a  function  of  the  incident  wave  height  wherever  the  cur¬ 
rent  is  strong.  The  effect  of  varying  the  incident  wave  height  on  H/H^ 
(defined  as  the  amplification  factor)  was  inves tiga ted ,  j n d  it  was  de¬ 
termined  that  the  greatest  variation  in  H/H^  occurred  in  the  deep  shadow 
zone  near  the  breakwater  and  shore  region  where  currents  were  strongest. 

In  both  Figures  22  and  23  (for  incident  wave  approach  directions  of  20 
and  30  deg,  respectively),  the  effect  of  initial  wave  height  is  seen  to 
be  greatest  near  the  breakwater  and  at  the  sensor  location  6  ft  from  the 
shoreline.  The  effect  of  incident  angle  is  also  apparent,  with  the  more 
acute  angle  producing  the  greater  degree  of  nonlinearity.  The  nonlinear 
effect  diminishes  rapidly  away  from  the  structure,  and  at  a  distance  of 
5  ft  downwave  appears  to  be  relatively  insignificant  except  for  the  ex¬ 
treme  shadow  region  at  acute  angles  of  incidence. 

Currents 


53.  The  establishment  of  a  counterclockwise  circulation  cell 
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a.  Distance  from  shoreline  =  6.00  ft  b.  Distance  from  shoreline  -  12.00  ft 

Figure  23.  Effect  of  varying  initial  wave  height  at  generator  on  wave-height  amplification,  H/Ho 
at  distances  of  6  and  12  ft  from  shoreline;  angle  of  incidence  -  30  deg 


approximately  4  ft  wide  adjacent  to  the  downwave  side  of  the  breakwater 
results  in  a  seaward  flowing  current  all  along  the  breakwater.  The 
bottom  current  is  especially  intense  approximately  40  percent  of  the 
breakwater  length  from  the  shoreline,  and  decreases  seaward  along  the 
structure  as  the  water  depth  increases  (since  mass  transport  remains 
constant)*  Middepth  velocities  decrease  seaward  along  the  breakwater 
also,  but  remain  approximately  35  percent  of  the  maximum  value  even  at 
the  seaward  end  of  the  structure.  The  middepth,  surface,  and  bottom 
velocity  measurements  are  displayed  in  Plates  49-51,  respectively.  The 
model  sidewall  boundary  is  responsible  for  the  creation  of  the  clockwise 
circulation  cell  that  develops  as  a  result  of  mass  transport  downcoast 
by  the  longshore  current.  The  numerical  model  of  Liu  and  Mei  (1975)  con¬ 
sidered  an  infinitely  long  coastline,  and  thus  precluded  the  clockwise 
circulation  pattern  which  developed  in  the  physical  model.  However,  the 
counterclockwise  circulation  expected,  and  indeed  which  occurred  adjacent 
to  the  model  breakwater,  d.d  not  appear  in  the  numerical  computat ional 
procedures.  A  counterclockwise  cell  developed  in  the  numerical  model  of 
Liu  and  Mei  (1975)  approximately  a  breakwater  length  downwave  from  the 
breakwater.  At  such  a  distance  in  the  physical  model,  however,  the 
sidewall  effects  begin  to  be  encountered.  Therefore,  the  physical  model 
results  for  a  finite  section  of  beach  should  not  be  expected  to  be 
identical  to  numerical  model  results  for  an  infinite  section  of  beach. 
Numerical  models  can,  however,  be  formulated  that  will  represent  the 
physical  model  boundary  conditions,  including  the  sidewall  effects.  The 
development  of  numerical  techniques  for  determining  wave  characteristics 
and  current  magnitudes  near  structures  is  presently  under  way.  These 
techniques  incorporate  combined  refraction  and  diffraction,  and  the  re¬ 
sults  of  these  experiments  will  be  used  to  verify  such  numerical  models. 
The  middepth  and  bottom  velocity  profiles  adjacent  to  the  breakwater  on 
the  downwave  side  are  shown  in  Figure  24  for  a  1-sec  wave  of  0.139-ft 
height  near  the  generator  and  approaching  from  an  angle  of  Incidence  of 
30  deg. 

54.  Coal  tracer  tests  were  conducted  to  obtain  qualitative 
knowledge  about  the  transport  and  dispersion  mechanisms.  A  1 /2-in .-thick 
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Middepth 
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of  incidence  *  30  deg 


layer  of  crushed  coal  was  placed  in  a  4-ft-wide  section  adjacent  downvave 
and  extending  the  length  of  the  breakwater.  The  coal  had  a  mean  diameter 
of  1.89  mm,  and  varied  in  size  from  a  maximum  diameter  of  3.35  mm  to  a 
minimum  diameter  of  0.43  mm.  The  model  was  operated  continuously  for 
1  hr  with  a  wave  period  of  1  sec,  and  the  results  are  displayed  in  Fig¬ 
ure  25.  Significant  movement  appeared  to  have  ceased  after  that  time, 
and  the  dynamic  balance  between  the  incoming  wave  train  and  the  return 
flow  created  by  the  circulation  cell  near  the  breakwater  resulted  in  the 
configuration  of  Figure  25b. 

55.  Net  fluid  particle  movement  was  visualized  by  a  series  of 
photographs  of  successive  locations  of  bottom  dye  streaks  originating 
from  four  locations  on  the  bottom  of  the  model.  The  dye  placed  at  a 
distance  of  4  ft  laterally  from  the  tip  of  the  breakwater  is  seen  in 
Photos  16-23  to  move  shoreward  and  become  part  of  the  development  of  the 
clockwise  circulation  cell  associated  with  the  downwave  boundary.  (A 
grid  spacing  of  2  ft  is  shown  in  these  photographs.)  Dye  placed  2  ft 
laterally  from  the  tip  of  the  breakwater  becomes  part  of  the  counter¬ 
clockwise  circulation  cell  near  the  breakwater.  After  the  circulation 
cells  had  developed  (Photos  21-23),  dispersion  of  the  dye  upward  from  the 
bottom  had  had  time  to  occur  and  the  streaklines  are  less  defined,  al¬ 
though  the  upward  dispersion  is  not  necessarily  associated  with  the  cell 
formation.  The  dye  movement  patterns  are  in  good  agreement  with  the 
bottom  velocity  measurements  shown  in  Plate  51.  A  circular  wave  pattern 
of  reflected  waves  radiates  outward  from  the  tip  of  the  breakwater,  and 
is  clearly  evident  in  this  series  of  photographs. 


PART  V:  SUMMARY  AND  CONCLUSIONS 


56.  The  purpose  of  this  study  was  to  obtain  detailed  and  precise 
experimental  data  regarding  wave-height  variations  and  currents  (patterns 
and  magnitudes)  downwave  from  a  shore-connected  breakwater  or  jetty 
under  the  simultaneous  effects  of  refraction  and  diffraction.  This  in¬ 
formation  provides  insight  into  the  phenomenon  of  combined  wave  refrac¬ 
tion  and  diffraction  and  can  be  used  to  verify  numerical  models  that 
simulate  this  phenomenon.  The  experimental  investigation  was  conducted 
in  a  wave  basin  that  was  molded  in  cement  mortar  and  consisted  of  an 
area  50  ft  x  60  ft  with  a  water  depth  of  1  ft  in  the  open-ocean  region. 

A  vertical,  impermeable  breakwater  (shore-connected)  extending  perpen¬ 
dicularly  from  the  shoreline  was  installed  on  a  beach  slope  of  IV  on  20H. 
Previous  experimental  studies  had  been  conducted,  for  the  most  part, 
with  horizontal  bottoms  which  provided  information  only  about  the  dif¬ 
fraction  phenomenon  (Blue  and  Johnson  1949,  Wiegel  1962).  One  experi¬ 
mental  study  of  extremely  limited  nature  (very  small  model,  72  sq  ft) 
had  been  conducted  that  propagated  waves  through  a  harbor  entrance 
(normal  angle  of  incidence)  onto  a  linearly  varying  bathymetry.  Results 
of  that  study  by  Mobarek  (1962)  indicated  that  as  long  as  the  incoming 
waves  were  within  the  transition-to-deepwater  range,  the  rule-of-thumb 
procedure  of  refracting  a  wave  to  the  structure,  diffracting  for  a  couple 
of  wavelengths  past  the  structure,  and  then  refracting  the  diffracted 
wave  shoreward  was  probably  adequate  for  most  applications.  However, 
for  long-period  (shallow  water)  waves,  the  combined  effects  of  refraction 
and  diffraction  should  be  considered  simultaneously. 

57.  This  experimental  study  used  34  wave-height  sensors  to  ascer¬ 
tain  the  average  maximum  wave  heights  downwave  of  the  shore-connected 
breakwater,  with  up  to  10  replications  of  each  individual  experiment 

to  quantify  the  experimental  variability.  The  measurements  were  highly 
accurate  with  a  standard  deviation  for  repeated  observations  of  usually 
2  to  3  percent  of  the  average  value  of  the  observations.  A  comparison 
of  these  experimental  data  of  combined  refraction  and  diffraction  was 
made  with  analytical  theories  of  diffraction  alone  for  a  shore-connected 


breakwater  on  a  horizontal  bottom.  As  expected,  the  comparison  was  found 
to  be  best  for  short-period  waves  and  for  waves  at  the  seaward  end  of 
the  breakwater,  since  short-period  waves  experience  less  shoaling  and 
refraction  and  these  effects  become  more  significant  the  farther  the 
waves  propagate  toward  the  shoreline.  The  incident  waves  generated  in 
the  model  were  sufficiently  small  that  the  waves  were  usually  linear 
within  the  measurement  region. 

58.  In  the  neighborhood  of  the  breakwater,  currents  existed  that 
affected  the  wave  heights.  The  magnitude  of  these  wave-induced  return 
currents  is  a  function  of  incident  wave  characteristics.  The  effect  of 
varying  the  incident  wave  height  on  the  wave-height  amplification  factor, 
H/Hq  ,  was  investigated  and  it  was  determined  that  the  greatest  variation 
in  H/Hq  occurred  in  the  deep  shadow  zone  near  the  breakwater  and  shore 
region  where  the  currents  are  the  strongest.  This  nonlinear  effect 
diminishes  rapidly  away  from  the  structure,  and  at  a  distance  of  5  ft 
downwave  appears  to  be  relatively  insignificant  except  for  the  extreme 
shadow  region  at  acute  angles  of  incidence.  The  establishment  of  a 
counterclockwise  circulation  cell  approximately  4  ft  wide  adjacent  to 
the  downwave  side  of  the  breakwater  results  in  a  seaward  flowing  current 
all  along  the  breakwater.  The  bottom  current  is  especially  intense  ap¬ 
proximately  40  percent  of  the  breakwater  length  from  the  shoreline  and 
decreases  seaward  along  the  structure.  The  physical  model  sidewall 
boundary  was  responsible  for  a  clockwise  circulation  cell  which  developed 
as  a  result  of  mass  transport  downcoast  by  the  longshore  current.  Based 
on  these  experiments,  any  numerical  techniques  for  describing  scour  and 
erosion  near  structures  should  account  for  wave-induced  currents  and 
circulation  cells  near  the  structure  as  well  as  tidal  or  other  currents 
that  may  exist  in  the  near  region. 
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Table  1 

Test  Conditions 


Incident  Wave  Approach  Direction 

30-deg 

Direction 

20- 

deg  Direction 

Wave  period. 

sec , 

« 

0.75 

1.00 

0.75 

1.00 

1.50 

Wavelength  near  generator. 

ft, 

2.82 

4.52 

2.82 

4.52 

7.69 

d/L 

m 

0.355 

0.221 

0.355 

0.221 

0.130 

Wave  height  near  generator. 

ft. 

- 

0.111 

0.122 

(Generator  stroke  *  0.50 

in. ) 

Wave  height  near  generator. 

ft. 

- 

0.164 

0.101 

0.180 

0.112 

(Generator  stroke  *  0.75 

in. ) 

Wave  height  near  generator. 

ft. 

- 

0.221 

0.139 

0.245 

0.154 

0.077 

(Generator  stroke  *  1.00 

in. ) 

Wave  height  near  generator. 

ft. 

- 

0.111 

(Generator  stroke  *  1.50 

in. ) 

Table  2 

Test  1 

Replications* 

Incident  Wave 

Approach  Direction 

30-deg 

20-deg 

Direction 

Direction 

Wave  period. 

sec , 

- 

0.75 

1.00  0 

.75 

1.00 

1.50 

Wavelength  near  generator. 

ft. 

- 

2.82 

4.52  2 

.82 

4.52 

7.69 

d/L 

* 

0.355 

0.221  0 

.355 

0.221 

0.130 

Wave  height  near  generator. 

ft. 

- 

0.111 

0 

.122 

Gage  arrangement 

1 ** 

10* 

Gage  arrangement 

2t 

10* 

5* 

Wave  height  near  generator. 

ft. 

- 

0.164 

0.101  0 

.180 

0.112 

Gage  arrangement 

1** 

10* 

10* 

3* 

3* 

Gage  arrangement 

2+ 

10* 

10* 

10* 

10* 

Wave  height  near  generator, 

ft. 

« 

0.221 

0.139  0 

.245 

0.154 

0.077 

Gage  arrangement 

!★* 

5* 

10* 

3* 

5* 

Gage  arrangement 

2t 

5* 

10* 

5* 

10* 

Wave  height  near  generator. 

ft. 

- 

0.111 

Gage  arrangement 

1  ** 

3* 

Gage  arrangement 

2+ 

*  Number  of  times  identical  tests  were  repeated  to  reduce  experimental 
variability. 

**  Gages  positioned  0,  2,  4,  6,  8,  10,  12,  and  14  ft  from  (normal  to) 
breakwater. 

t  Gages  positioned  1,  3,  5,  7,  9,  11,  13*  and  15  ft  from  (normal  to) 
breakwater. 


Table  3 


Distribution  of  Wave  Energy  Throughout  the  Model,  30-deg  Incident 


Direction,  Period  =  0.75  sec.  Ocean  Wave  Height  =  0.111  ft 


Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Fundamental 

First 

Period 

Harmonic 

Gage 

0.75  sec 

0.38  sec 

1* 

99.91(-162) 

0.08 (+25) 

2* 

99. 8l(— 177) 

0. 18(+27) 

3 

99.290151) 

0. 70(+60) 

4 

98 . 74 (+74) 

1 . 23(+149) 

3 

99 . 44 (—57) 

0. 55(-108) 

6 

99 . 00(-170) 

0. 99(+23) 

7 

99. 40067) 

0 . 60(+140) 

8 

99. 2  6 (— 6  9 ) 

0. 73(-150) 

9 

99.58(+159) 

0. 42(-19) 

10 

99.96076) 

0. 03(-141) 

11 

99.23094) 

0. 76C-171) 

12 

98 . 7  5 (-35) 

1 . 24 (-67) 

13 

99.37 (-165) 

0. 62(+41) 

14 

98 . 50086) 

1.46(+161) 

15 

98.25(-60) 

1 . 68(-112) 

16 

98 . 95 (+168) 

1.03(-32) 

17 

99.34  065) 

0. 66(+137) 

18 

99 . 9602) 

0.0H+33) 

19 

97 . 15012) 

2. 77 (-25) 

20 

97.690128) 

2. 24 (+100) 

21 

96.860115) 

3 . 03(-134) 

22 

97 .41016) 

2. 54(-38) 

23 

97.950156) 

1. 99(+44) 

24 

99 . 05084) 

0. 94(+149) 

25 

98.87 (-9) 

1 . 13 (-24) 

26 

99 . 84 (-59) 

0. 16(-152) 

27 

94 . 67  068) 

4 . 90(-146) 

28 

92.340175) 

6. 96(-l) 

29 

91 . 37 (+45) 

7 . 77 (+76) 

30 

92.16088) 

7 . 22 (+1 75) 

31 

93 . 920136) 

5 . 67 (-96) 

32 

98 . 10037) 

1 . 87 (+73) 

33 

98. 61054) 

1 . 37 (-120) 

34 

99.  78090) 

0.21(+172) 

Second  Third  Fourth 

Harmonic  Harmonic  Harmonic 

0.25  sec  0.19  sec  0.15  sec 

0.01 (+19)  O.OOC-137)  0.00(-170) 

0. 01 (+12)  0. 00(-136)  0. 00(+59) 

0. 01 (-87)  0. 00 (+103)  0. 00(-64) 

0. 02(-112)  0. 00(-9)  0.00(+101) 

0.01(-102)  0. 00(-167)  0. 00(-32) 

0.0K-116)  0. 00(+61)  0. 00(+167) 

0. 00(-53)  0. 00(+147)  0.00(-50) 

0.01(+144)  0. 00 (+8 3)  0. 00(+73) 

0. 00(-158)  0. 00 (-154)  0.00(-44) 

0.01(+17)  0. 00(-149)  0. 00(-45) 

0.01(-59)  0. 00(+40)  0. 00(+121) 

0. 12(-64)  0. 00(-142)  0.00(+154) 

0.01(-95)  0. 00 (+3)  0. 00(-178) 

0. 03(-115)  O.Ol(-ll)  0. 00 (+8 3) 

0. 06(-152)  0.01(+170)  0. 00(+165) 

0.01(+130)  0. 00(-64)  0.00(+110) 

0. 00 (-23)  0. 00(+113)  0. 00(-154) 

0. 02(+108)  0.00(-58)  0.00(-76) 

0. 08 (-35)  0. 00(-28)  0.00(+75) 

0. 07 (-38)  0. 0H-149)  0.00(+152) 

0. 10(-17)  0 . 00 (+118)  0. 00(-40) 

0. 05 (-52)  0. 00(-48)  0.00(+25) 

0. 05 (-109)  0. 00(+119)  0.00(+ll) 

0. 01 (-176)  0. 00(-119)  0.00(-5) 

0. 00(-48)  0. 00(-149)  0.00(-148) 

0. 00(+99)  0. 00(+23)  0.00(-63) 

0. 36(+140)  0. 06 (+8 7)  0.01(+61) 

0. 62 (+178)  0. 08 (+21)  0.00(-115) 

0. 77 (+119)  0. 09(+176)  0.00(-129) 

0. 57 (+85)  0.04 (+] 3)  0.00(+25) 

0. 37 (+38)  0. 04(-174)  0.00(+12) 

0. 03 (+126)  0. 00(+113)  0. 00(+119) 

0. 0H-145)  0. 00(+93)  0.00(-0) 

0. 0H-144)  0. 00(-130)  0. 00(-145) 


*  Ocean  gage. 


r 


Table  4 

Distribution  of  Wave  Energy  Throughout  the  Model,  30-deg 

Direction,  Period  =  0.75  see.  Ocean  Wave  Height  -  0.164  ft 
Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Gage 

Fundamental 

Period 

0.75  sec 

First 
Harmonic 
0.38  sec 

1* 

99.62(+67) 

0. 30(+158) 

2* 

99. 93 (+58) 

0.03(+31) 

3 

98. 22(-29) 

1.68(-46) 

4 

98 . 64 (-143) 

1 . 32(+58> 

5 

97 . 41(+88) 

2 . 40(-170) 

6 

96 . 87 (-12) 

2 . 91 (-49) 

7 

96. 18(-132) 

3 . 39(+112) 

8 

97 . 62 (+109) 

2. 15(-167) 

9 

99. 20(+4) 

0. 7  3 (—40) 

10 

99.03(-52) 

0. 78(-35) 

11 

96. 10(-131) 

3 . 66(+92) 

12 

97. 84 (+108) 

2 . 10(-141) 

13 

94 . 61(-7) 

4 . 92(-19) 

14 

97 . 26(-121) 

2. 63(+114) 

15 

94 . 56(+116) 

4 . 99(-130) 

16 

98 . 09(+5) 

1. 79(+3) 

17 

98 . 7  2 (—84) 

1. 16(+144) 

18 

99. 96(-117) 

0. 04(-32) 

19 

96.06(+133) 

3 . 69(-107) 

20 

92. 96(+22) 

6. 43(+39) 

21 

94 . 45(-103) 

5. 28 (+146) 

22 

90. 03(+132) 

8 . 89(-100) 

23 

95. 28(+14) 

4.42(+18) 

24 

97 . 22(-81) 

2. 63(-178) 

25 

98. 02(-152) 

1.94 (+14) 

26 

99. 64(-169) 

0. 31(-54) 

27 

87. 58 (+65) 

10. 83(+126) 

28 

86. 90(-51) 

11. 26(-114) 

29 

80. 76(+170) 

15. 7  3 (+5  6) 

30 

85. 36(+60) 

12.47(+106) 

31 

89.89(-49) 

9. 07 (-102) 

32 

95.04(-U7) 

4. 67 (+101) 

33 

97.63(+173) 

2. 25 (—63) 

34 

98. 69(+171) 

1. 30(-59) 

Second 
Harmonic 
0.25  sec 

Third 
Harmonic 
0.19  sec 

Fourth 
Harmonic 
0.15  sec 

0. 03(+23) 

0. 02(+128) 

0 . 02 (-156) 

0. 00 (+116) 

0. 03 (-118) 

0. 00(-49) 

0.09(-52) 

0. 01 (-41) 

0.00 (+7) 

0. 03 (-98) 

0. 01 (-174) 

0. 00(+170) 

0. 17 (-56) 

0.02 (+85) 

0. 00 (-125) 

0. 20(-81) 

0. 02(-lll) 

0. 00(+169) 

0. 36 (+17; 

0. 07 (-87) 

0. 00 (+175) 

0. 21 (-72) 

0.02 (+62) 

0. 00(+74) 

0. 06 (-133) 

0. 02(-173) 

0.00(-174) 

0. 17 (-48) 

0. 02(-79) 

0. 00 (-95) 

0. 20(-38) 

0. 04(-149) 

0. 00(+99) 

0.05 (+4) 

0. 01(+166) 

0. 00(-24) 

0. 43 (-22) 

0. 04(-15) 

0. 00 (-71) 

0. 09(+13) 

0. 01(-60) 

0. 00 (-36) 

0. 40(-4) 

0. 05(+139) 

0.00 (+68) 

0.11(-1) 

0.01(-4) 

0.00(+101) 

0. 10(+5) 

0. 02(-92) 

0. 00(+161) 

0. 00(+45) 

0. 00(-121) 

0. 00(+43) 

0.  23 (+20) 

0.  02 (-165) 

0. 00(+113) 

0.57 (+74) 

0.04 (+116) 

0. 00 (+112) 

0. 26(+42) 

0. 01 (-24) 

0.00(+lll) 

1. 00(+32) 

0.07 (+172) 

0.  01 (-98) 

0.27 (+30) 

0. 02(+65) 

0. 00(-27) 

0.14 (+87) 

0.01 (+10) 

0. 00(-78) 

0. 04 (-175) 

0. 00(+7) 

0.  00(-171) 

0.04 (+81) 

0.01(-30) 

0.00 (+121) 

1. 46(-171) 

0. 14 (-98) 

0. 00(-46) 

1. 65(-171) 

0.16 (+128) 

0.04 (+35) 

2. 86 (+134) 

0.  50 (—67) 

0. 15 (+7  7 ) 

1.95 (+164) 

0. 19(-137) 

0. 03 (-114) 

0.91(-149) 

0. 13(— 177) 

0. 00 (+17 6) 

0. 26(-40) 

0. 03(-173) 

0. 00(+77) 

0.11 (+70) 

0. 01 (+178) 

0.00(-l) 

0.01 (+31) 

0.01 (+129) 

0.  00(-33) 

*  Ocean  gage. 


Table  5 


1) i s t r ibution  of  Wave  Energy  Throughout  the  Model,  30-deg  Incident 
Direction,  Period  =  Q.75  sec,  Ocean  Wave  Height  =  0* 221  ft 
Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Fundamental 

F  irs  t 

Period 

Ha rmonic 

Gage 

0.75  sec 

0.38  sec 

1* 

98 . 03 (-87) 

1 . 63(-141) 

2* 

99.49(-73) 

0.48(-179) 

3 

96.62(+l?7) 

3 . 25 (-12) 

4 

92 . 90 (+41) 

6.44(+77) 

5 

95 . 82(-98) 

3.  79 (  —  17  3) 

6 

94 . 56(+151) 

4 . 92(-73) 

7 

94 . 24 (+20) 

5.31(+52) 

8 

95.02(-102) 

4 . 60(+148) 

9 

96 . 56(+164) 

3 .18(-55) 

10 

97. 73(+118) 

2. 19(-112) 

11 

92.23(+55) 

6 . 80(+113) 

12 

96.05(-76) 

3 . 81(-146) 

13 

91.37(+152) 

7 . 61(-65) 

14 

91.71(4-27) 

7.47 (+56) 

15 

92.90(-101) 

6. 57 (+149) 

16 

94.90(4-154) 

4. 62(-70) 

17 

97.17(4-75) 

2 . 66(+137) 

18 

96.88(4-48) 

2.86(+117) 

19 

93 . 00(-54) 

6 . 48(-106) 

20 

87.63(4-174) 

10. 82(-27) 

21 

85.85(4-31) 

11. 86(+62) 

22 

89.21 (-88) 

9. 63(+170) 

23 

92. 98 (+149) 

6.48(-62) 

24 

95 . 48 (+74) 

4 . 1 9 ( +1 22) 

25 

96. 3  2 (+2 ) 

3 . 50(-ll) 

26 

95 . 7  9  ( —  3) 

3.86(+13) 

27 

83. 14(-145) 

14. 17 (+60) 

28 

78. 60(+81) 

16. 72(+l42) 

29 

82.84(-62) 

1 4 . 74(-145) 

30 

77.44(-168) 

17 . 58 (+9) 

31 

84. 4 5 (+8 2) 

13.40(+154) 

32 

94 . 29(+26) 

5 . 32(+54) 

33 

94 . 35 (-26) 

5. 20(-68) 

34 

92. ll(-29) 

7 . 36(-60) 

Second 
Harmonic 
0.25  sec 

Third 
Harmonic 
0.19  sec 

Fourth 
Harmonic 
0.15  sec 

0. 29 (-174) 

0. 03(+32) 

0. 12(+172) 

0. 60(+123) 

0. 35 (+136) 

0. 05 (+105) 

0. 00(+103) 
0.01 (+58) 

0. 05 (+165) 
0.04 (+80) 

0. 00(+105) 
0. 00(-39) 

0. 00(+48) 
0.01(+163) 
0. 00(+44) 

0.47 (+65) 

0. 40(+105) 

0. 35(+53) 

0. 23(+108) 

0. 06 ( —  5) 

0.04(-155) 

0. 05(+157) 

0. 01 (-46) 
0.03(-53) 
0.01 (+128) 

0 . 01 (-43) 
0.00 (+122) 
0.01 (+150) 
0.00(-154) 
0. 00(-72) 

0.84 (+175) 

0. 1U-178) 
0.87 (+83) 

0. 72(+103) 

0. 48 (+4 7) 

0. 11(-119) 

0. 03 (+1 73) 

0. 11  (-137 ) 

0. 08(+143) 

0. 03(-7 ] ) 

0.03(-92) 
0.00(-100) 
0.03(-3) 
0.02 (+128) 
0.02 (+153) 

0.  45 (+78) 

0. 13 (-157) 

0. 24(+179) 

0. 49 (-148) 
1.33 (+137) 

0. 03(-120) 

0. 04 (-39) 

0. 02 (-120) 

0. 03 (-166) 

0. 17 (-72) 

0. 00(-41) 
0.00(+54) 
0. 01 (-60) 
0.  00(+l 5) 
0. 06(+84) 

1 . 80(+94) 

1. 06 (+74) 

0. 49(+102) 

0.  31 (+171) 

0. 16(-23) 

0.04(+ll7) 
0.08 (-32) 

0. 04 (-99) 

0. 02 (-84) 

0.  01 (-14) 

0.01 (+137 ) 
0.02 (-159) 
0.01 (+12) 
0. 00(+7) 
0.00(+15) 

0. 29(+34) 

2.  26(-98) 

3. 57 (-159) 

1 . 98(+125) 

3. 82(-171) 

0.05(+46) 

0. 34 (+94) 

0.  79(-lll) 
0. 33(+15) 

0. 88 (+2) 

0.01 (+51) 
0. 09(-8] ) 
0. 32(-67) 
0.11 (-91) 
0.  28(+167) 

1 . 87 (-123) 
0.  36(+99) 

0. 40(-109) 

0. 44 (-84) 

0.  22(-44) 

0. 03(+156) 
0.05(-127) 
0.08 (-96) 

0. 05 (+8) 

0. 00(-91 
0.01 (-125) 
0. 01 (-95) 

*  Ocean  gage. 


Table  6 


Distribution  of  Wave  Energy  Throughout  the  Model,  30-deg  Incident 
Direction,  Period  =  1.00  sec.  Ocean  Wave  Height  =  0.101  ft 
Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Gage 

Fundamental 

Period 

1 , 00  sec 

First 
Harmonic 
0.50  sec 

Second 
Harmonic 
0.33  sec 

Third 
Harmonic 
0.25  sec 

Fourth 
Harmonic 
0. 20  sec 

1* 

2* 

3 

4 

5 

98.87 (+130) 
99.42(+134) 
98. 96(+129) 
97. 95 (+49) 
97.52(-31) 

0. 77 (+178) 

0. 23 (+17  4) 

1 . 03 (-157) 

2 . 02(+127) 

2 . 46(-99) 

0. 23 (-125) 

0.  22 (-145) 

0.  00(-126) 

0. 02(+167) 
0.02(-178) 

0. 11(-14) 

0. 11(-18) 
0.00(-42) 

0. 00(+95) 

0. 00(+78) 

0.02 (+114) 
0.03 (+105) 
0. 00(+86) 

0 . 00 (+178) 
0.00 (+9) 

6 

99 . 66(-lll) 

0. 26(+147) 

0.07 (+159) 

0. 01(+58) 

0. 00(+2) 

7 

97. 07 (+172) 

2 . 84 (-5) 

0. 09 (-155) 

0.01 (+66) 

0. 00 (-144) 

8 

97 . 96(+88) 

1 . 95 (+160) 

0. 06 (+144) 

0. 02(-164) 

0. 01 (-85) 

9 

99.43(+27) 

0.43(-13) 

0. 13(+14) 

0.00 (+44) 

0. 00(-84) 

10 

99. 9  0 ( +4 ) 

0. 10(+115) 

0.00(+110) 

0. 00(+157) 

0. 00(-178) 

11 

95 . 53 (-57) 

4.41(-133) 

0. 05 (+170) 

0.01(+95) 

0. 00(+61) 

12 

99. 61(-139) 

0 . 37 (+54) 

0. 01 (-159) 

0. 00(+52) 

0. 00 (-21) 

13 

96. 30(+140) 

3 . 59(— 7  6) 

0. ll(+49) 

0. 00(+169) 

0. 00 (-56) 

14 

96.58(+60) 

3 . 26(+89) 

0 . 12 (+109) 

0 . 03 (+133) 

0.00(-178) 

15 

98. 95(-19) 

1 . 03 (-46) 

0. 01 (-161) 

0.01(+172) 

0. 00(+125) 

16 

99.32(— 97) 

0 . 61 (-174) 

0. 06(-169) 

0.01(+119) 

0. 00(+65) 

17 

97 . 47 (-147) 

2. 49(+78) 

0.01(+107) 

0. 03(~28) 

0. 00(-135) 

18 

98. 12(-160) 

1 . 85 (+12) 

0. 02(+167) 

0. 00(-134) 

0. 00(+69) 

19 

95.16(+121) 

4 . 65 (-111) 

0. 17 (+27) 

0 . 01 (+150) 

0. 00 (-7 9) 

20 

95. 02(+51) 

4 . 66(+70) 

0.  29 (+111) 

0.  03 (+154) 

0. 00 (-144) 

21 

96 . 02(-36) 

3.82(-80) 

0. 16(-126) 

0. 01 (+176) 

0.  00(-160) 

22 

92. 75(-116) 

6 . 70(+117) 

0.  50(-7) 

0.04 (-142) 

0.01(+119) 

23 

95 . 80(+164) 

4.01(-52) 

0. 17 (+111) 

0. 01 (-108) 

0.01 (+69) 

24 

97.66(+90) 

2 . 03 (+145) 

0. 30(-130) 

0.01(-44) 

0. 00(+20) 

25 

99 . 64 (+5  5) 

0 . 33(+96) 

0. 02 (—  39) 

0.01(+4) 

0. 00(+115) 

26 

99 . 67 (+45) 

0.32(+114) 

0.01 (+20) 

0. 01 (+31) 

0.00(+110) 

27 

89 . 24 (-14) 

9 . 52 (-51) 

1. 08(-81) 

0. 14 (-97) 

0. 03 (—97 ) 

28 

82. 63(-96) 

14 . 39(+142) 

2. 48 (+33) 

0. 43 (-74) 

0.  07 (-171) 

29 

87.19(+174) 

11.1K+38) 

1.41 (+118) 

0. 25 (-78) 

0.04 (+96) 

30 

88 . 20(+92) 

10. 37 (+168) 

1. 21 (-109) 

0. 18 (-23) 

0.04(+83'> 

31 

90.84 (+16) 

8 . 36(+24) 

0. 72(+27) 

0 . 07 (+31) 

0.01 (+47) 

32 

91 . 39 (— 32) 

8 . 40(-78) 

0. 20(-138) 

0.  00(+l42) 

0.  01 (+85) 

33 

93.99(-74) 

5 . 79 (-173) 

0. 19 (+98) 

0.  02 (-35) 

0. 01 (-179) 

34 

96. 13(-85 

3.75(+157) 

0. 1 2 (+64 ) 

O.OO(-IOO) 

0.00 (+96) 

*  Ocean  gage. 


Table  7 

Distribution  of  Wave  Energy  Throughout  the  Model,  30-deg  Incident 
Direction,  Period  =  1.00  sec,  Ocean  Wave  Height  =  0.139  ft 
Percent  of  Total  Energy  and  Phase  Angle  (deg) 


(Sage 

Fundamental 

Period 

1,00  sec 

First 

Ha rmonic 
0.50  sec 

Second 
Harmonic 
0.33  sec 

Third 

Harmonic 

0.  25  sec 

Fourth 
Harmonic 
0.20  sec 

1* 

99. 23(+127) 

0. 65(+126) 

0. 09(+124) 

0.  03(-115) 

0. 01 (+52) 

2* 

99. 77 (+126) 

0.14 (+31) 

0.  05 (+2 2) 

0.  03 (+141) 

0. 00(-41) 

3 

97. 70(+lll) 

2. 27 (-100) 

0.03 (+15) 

0. 00 (-115) 

0. 00(-139) 

4 

94. 44 (+35) 

5. 30(+53) 

0. 25 (+60) 

0.02 (+76) 

0.00 (+134) 

5 

99. 25 (-43) 

0. 64(-lll) 

0. 08 (+63) 

0.  03(+ll) 

0.00 (+7) 

6 

96. 35 (-120) 

3. 55(+133) 

0.08 (+61) 

0.01 (-4) 

0. 00(-86) 

7 

96. 15(+164) 

3. 77  (-51) 

0.08 (+108) 

0. 00(-52) 

0. 00(+121) 

8 

99. 53(+76) 

0. 36 (+123 ) 

0.09 (+43) 

0.  02 (+140) 

0. 00(-148) 

9 

98. 44 (+10) 

1 . 53 (+38) 

0. 02 (-154) 

0.01 (-104) 

0. 00(-56) 

10 

99. 39 (— 14 ) 

0.58 (+3) 

0.03 (-53) 

0.  00 (-25) 

0.00(-86) 

11 

98.51(-75) 

1.40(+172) 

0.08 (+76) 

0.01 (+29) 

O.OO(-ll) 

12 

96. 78 (-157) 

3. 02(+68) 

0.19(-84) 

0.01 (+160) 

0. 00 (+7 ) 

13 

92 . 87 (+130) 

6. 39(-121) 

0. 66 (-2) 

0. 08(+122) 

O.Ol(-lOO) 

14 

97 . 63 (+49) 

2. 27 (+84) 

0.09 (+98) 

0.01 (+146) 

0. 00 (-86) 

15 

95. 76(-32) 

4.03(-46) 

0.  20 (-82) 

0. 02 (-129) 

0. 00(-89) 

16 

96. 65 (—111) 

3. 08 (+125) 

0. 23 (+28) 

0.04 (-75) 

0. 01 (-162) 

17 

98. 53(-168) 

1 .42 (+28) 

0. 04 (-93) 

0.01 (+156) 

0. 00(+70) 

18 

99.06(+178) 

0.92 (-15) 

0. 02 (+129) 

0.  00(-113) 

0. 00(+105) 

19 

90. 76(+104) 

8.37 (-178) 

0. 73 (-77) 

0.11 (+37 ) 

0.02 (+163) 

20 

95. 46(+31) 

4. 18 (+48) 

0.  33 (+7 7) 

0. 04 (+116) 

0. 00 (+176) 

21 

88. 63(-51) 

10.09(-112) 

1. 13(-167) 

0. 14 (+145) 

0.01 (+103) 

22 

87 . 97 (-127) 

10. 32 (+79) 

1.42(-58) 

0.  25 (+168) 

0.03 (+38) 

23 

96 . 86 (+150) 

3. 06(-76) 

0. 06 (+79) 

0.01(-109) 

0. 01 (+83) 

24 

96 . 65 (+76) 

3. 34 (+150) 

0. 01 (-94) 

0. 00 (+] 59) 

0. 00(-134) 

25 

97. 80(+37) 

2 . 13 (+71) 

0.07 (+86) 

0. 00  (+4  5) 

0.00 (+82) 

26 

97.91 (+24) 

2. 08 (+57) 

0.01 (+9) 

0.  00(+115) 

0.  00 (-43) 

27 

73. 48 (— 39 ) 

20. 34 (-108) 

4. 83 (-161 ) 

1 . 14 (+148) 

0.21 (+95) 

28 

78. 23 (—121 ) 

17.16 (+88) 

3. 85(-52) 

0.  68(+171) 

0. 08 (+21) 

29 

76. 3  5 (+1 51 ) 

18. 64 (+112) 

4.07 (+70) 

0. 81 (-142) 

0. 1 3 (— 1 ) 

30 

74 . 1 3 (+7  4 ) 

20. 00(+130) 

4. 57 (-165) 

10. 96(-97) 

0. 21(-29) 

31 

83. 78 (+4 ) 

14 . 30(-19) 

1. 65(-32) 

0.24 (-38) 

0. 03 (-32) 

32 

92. 30(+49) 

6. 95 (— 121) 

0. 65(+164) 

0.  09(+103) 

0.01 (+45) 

33 

94. 90(-94) 

4 . 92 (+145) 

0. 17 (+30) 

0.  00(-63) 

0. 00(+135) 

34 

95. 88(-103) 

3. 89(+126) 

0. 23 (-19) 

0.01 (-124) 

0. 00 (+31) 

*  Ocean  gage. 
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Table  8 

Distribution  of  Wave  Energy  Throughout  the  Model,  20-deg  Incident 
Direction,  Period  =  0.75  sec.  Ocean  Wave  Height  =  0.122  ft 
Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Gage 

Fundamental 

Period 

0.73  sec 

First 
Harmonic 
0.38  sec 

Second 
Harmonic 
0.25  sec 

Third 
Harmonic 
0.19  sec 

Fourth 

Harmonic 

0.15  sec 

1* 

99 . 93 (+110) 

0.02(4-102) 

0.02(4-61) 

0.02(4-102) 

0 . 00 (-180) 

2* 

99 . 57 (+149) 

0. 36 (-108) 

0.05(4-23) 

0. 02 (-170) 

0. 00(-91 ) 

3 

98 . 14 (-119) 

1.74(4-163) 

0.10(4-86) 

0.0H-2) 

0.00(-79) 

4 

99.09(4-161) 

0 . 86 (-61 ) 

0.04(4-64) 

0.01 (-118) 

0.00(475) 

3 

97.86(4-67) 

2.07(4-150) 

0. 07 (-113) 

0.01 (-0) 

0.01(4131) 

6 

98.42 (-21) 

1 . 54 (-59) 

0.03(-114) 

0.00(-144) 

0.00(477) 

7 

99.63 (-101) 

0.37(4-147) 

0.00(-105) 

0.00(4-122) 

0.00 (+48) 

8 

97.26(4-164) 

2 . 64 (-24 ) 

0.10(4-148) 

0. 00(-28) 

0.00(-108) 

9 

98.49(4-70) 

1.34  (4-118) 

0.14(4-172) 

0.03(-103) 

0.00(-l) 

10 

99.88(4-10) 

0.09(4-105) 

0. 02 (-135) 

0 . 00(-58) 

0.00(-63) 

11 

97.81(4-157) 

2. 05(-38) 

0.13(4-148) 

0.01 (-49) 

0.00(4126) 

12 

98.83(4-67) 

1.16(4-114) 

0.01(4-125) 

0.00(-l 59) 

0.00 (-105) 

13 

97 . 64 (-18) 

2 . 27 (-26) 

0. 09(-28) 

0.01 (-28) 

0.00(430) 

14 

97 . 77 (-106) 

2.11(4-129) 

0.10(4-1) 

0.01 (-93) 

0.00(4179) 

13 

99.24(4-170) 

0. 76(-27) 

0.01(4-123) 

0. 00(-22) 

0.00(4178) 

16 

97.11(4-66) 

2.70(4-145) 

0. 1 7 (-139) 

0.03(-64) 

0.00(491) 

17 

97. 31 (-15) 

2 . 64 (-43) 

0.05 (-65) 

0. 00(-60) 

0.00(-67) 

18 

99 . 69 (-60) 

0.29 (-171) 

0.02(4-24) 

0.01 (-49) 

0.00 (-81) 

19 

95.27(4-73) 

4.43(4-1  38) 

0.27 (-148) 

0.02 (-41) 

0. 00 (-76) 

20 

97 . 95 (-11 ) 

2.01 (-25) 

0 . 03 (-14) 

0.01(-47) 

0.00(481) 

21 

96. 20 (—  99 ) 

3.  58(4-162) 

0.20(4-72) 

0.01 (-9) 

0.00 (-7 5) 

22 

96.70(4-174) 

3 . 19 (-28) 

0.10(4-132) 

0.01 (-54) 

0.00(4162) 

23 

97.67(4-83) 

2.25(4-166) 

0.08 (-100) 

0.00(4-26) 

0.00(4113) 

24 

96 . 28 (-22) 

3. 49 (-43) 

0. 20 (—6 1 ) 

0.03(-64) 

0.00(-22) 

23 

97 . 24 (-86) 

2.64(4-175) 

0.12(4-72) 

cr- 

1 

o 

o 

0.00(-102) 

26 

98. 96(-119) 

1.03(4-80) 

0.01(4-86) 

0. ()()(-!  51 ) 

0.00(480) 

27 

91.29(4-27) 

8.03(4-35) 

0.62(4-54) 

0.05(4-92) 

0.00(4179) 

28 

91 . 63 (-57) 

7.45(— 118) 

0.80(-172) 

0.11(4-144) 

0.00(4115) 

29 

90.62(4-147) 

8.43(4-89) 

0.75 (-92) 

0.09(4-133) 

0.11(41) 

30 

92.  75(4-123) 

6 . 83 (-120) 

0.39(4-2) 

0.03(4-140) 

0.00(-93) 

31 

89.96(4-23) 

8.97(4-41) 

0.94  (4-67) 

0.12(4-102) 

0.00(4143) 

32 

93 . 64 (-66) 

5 . 8  3  ( — 14 1 ) 

0.48(4-147) 

0.05(493) 

0.00(465) 

33 

96. 12 (-122) 

3.68(4-88) 

0. 19 (-50) 

0.01 (4171) 

0.00(458) 

34 

98.42(-147) 

1.56(4-37) 

0.02(-157) 

0.00(459) 

0.00 (-35) 

*  Ocean  gage. 
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Table  9 

Pis t ribut i o n  of  Wave  Energy  Throughout  the  Model,  20-deg  Incident 
Direction,  Period  =  0.75  sec,  Ocean  Wave  Height  =  0.180  ft 
Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Gage 

Fundamental 

Period 

0. 75  sec 

First 
Harmonic 
0.38  sec 

1* 

99. 27 (+68) 

0.63 (+70) 

2* 

96. 88(+107) 

2. 90(-171) 

3 

98. 06(+124) 

1 . 83(-147) 

4 

96. 37 (+21) 

3. 22 (+50) 

5 

97.47 (-53) 

2. 37 (-149) 

6 

99. 19(-143) 

0.68 (+97) 

7 

95. 15(+139) 

4. 34 (-71) 

8 

94.67 (+56) 

4 . 78 (+93) 

9 

97. 00 (-17) 

2. 60(-86) 

10 

99.48 (-41) 

0. 50(-125) 

11 

97.77 (+14) 

2. 09 (+44) 

12 

93. 88 (-68) 

5. 56(-147) 

13 

98. 40(-151) 

1.57 (+40) 

14 

94. 56(+123) 

4 . 95 (-99) 

15 

92. 82 (+42) 

6.44 (+73) 

16 

96. 17 (-41) 

3. 52 (-100) 

17 

96. 95 (-100) 

2. 75(+108) 

18 

98. 56 (-112) 

1.39 (+106) 

19 

91. 23(-60) 

7.86(-125) 

20 

95. 25  (-154) 

4. 59 (+4 7) 

21 

93. 72 (+121) 

5. 70(-lll) 

22 

90. 49 (+40) 

8. 23 (+62) 

23 

94.48(-45) 

5. 16(-113) 

24 

95. 12(-125) 

4.54 (+86) 

25 

97. 77 (-170) 

2. 09(-21) 

26 

97. 51 (-168) 

2. 32 (+12) 

27 

88.47 (-120) 

10. 33(+105) 

28 

82. 35(+145) 

14. 75(-78) 

29 

84. 25 (+65) 

13. 37 (+109) 

30 

84.42(-35) 

13. 26 (—  75 ) 

31 

85. 2  3 (—114 ) 

12. 69(+l 30) 

32 

94. 57 (-173) 

5. 08 (+6) 

33 

97 . 22 (+159) 

2. 65(-78) 

34 

95. 33(+160) 

4.46(-57) 

Second 

Harmonic 

0. 25  sec 

Third 

Harmonic 

0.  19  sec 

Fourth 
Harmonic 
0. 15  sec 

0.04 (+59) 

0. 18 (-54 ) 

0. ll(-66) 

0.  34 (+101) 

0. 16 (+131) 

0. 05(-3) 
0.03 (+99) 

0. 00 (+41) 

0. 06(+135) 

0. 00(+46) 

0.01 (+49) 

0. 01 (-129) 
0. 00(+52) 

0.00(+11m) 

0, 00(-77) 

0. 09(+49) 

0. 45 (+86) 

0. 50(+144) 

0. 36 (-125) 
0.02 (+124) 

0. 02 (-68) 

0. 06(-112) 

0. 05(-137) 

0. 04 (-126) 
0.01 (-145) 

0.02(-147) 
0. 00(+105) 
0. 00(-136) 
0. 00 (-120) 
0. 00 (+155) 

0. 12 (+103) 

0. 52 (+14 1 ) 

0. 02(-168) 
0.41 (+59) 

0. 67 (+119) 

0. 02 (+161) 

0. 05 (+86) 

0. 01 (-115) 

0. 06(-147) 

0. 07 (+172) 

0. GO(-87) 

0. 00(-29) 

0. 00 (+158 
0.00 (+63) 
0.01 (+174) 

0. 29 (-147) 

0.  26  (-19) 

0. 05(-64) 

0. 81 (-177) 
0.16 (-108) 

0. 02 (-173) 

0. 04(-95) 
0.01 (+150) 

0. 09(+145) 

0. 00(+157) 

0.00 (+63) 

0. 00(-176) 
0.  00 (+9) 

0. 00(+60) 

0. 00 (+176) 

0. 52 (+30) 
1.10 (+96) 
0.34 (+179) 

0. 33 (—58 ) 

0. 14 (+124) 

0. 06 (+170) 

0. 14 (+128) 

0. 01 (+105) 
0.01 (-161) 

0. 01 (-37) 

0. 00 (-107 ) 
0.03(+133) 
0.01 (-6) 

0. 00(-87) 

0. 00(-172) 

0. 13(-178) 

1. 10(-20) 

2. 49 (+64) 

2. 09 (+159) 

1 . 98  (-111) 

0. 03(+14) 

0. 09(-136) 

0. 34 (-159) 

0. 23  (—161) 
0.27 (-149) 

0. 01 (-139) 
0.01 (+60) 

0. 08(-37) 

0. 06 (-139) 
0.07(+154) 

1. 88 (+23) 

0. 31 (-163) 

0. 1 2 (+31 ) 
0.18 (+84) 

0. 19 (-79) 

0. 03 (+75) 

0. 00(+175) 

0. 02(-119) 

0. 02(+138) 
0.01 (-33) 

0. 00(+26) 

0. 01 (+89) 

*  Ocean  gage. 


Table  10 

Distribution  of  Wave  Energy  Throughout  the  Model,  20-deg  Incident 
Direction,  Period  =  0.75  sec.  Ocean  Wave  Height  =  0.245  ft 


Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Gage 

Fundamental 

Period 

0.75  sec 

First 
Harmonic 
0.38  sec 

1* 

97 . 71 (-79) 

1 . 93 (-103) 

2* 

95.44(-24) 

3. 74 (-79) 

3 

98. 59 (+19) 

1 . 40 (+12) 

4 

90. 28 (-89) 

8. 62 (+1 77 ) 

3 

97. 92(+170) 

1 . 99 (-9) 

6 

88. 98 (+83) 

9. 50 (+165) 

7 

94. 42 (+2) 

5. 10(-18) 

8 

94. 81 (-84) 

4.84 (-149) 

9 

94. 37 (-167) 

5. 05 (+4 9) 

10 

95. 39(+139) 

3. 75(-44) 

11 

89. 71 (-92) 

8. 89 (+175) 

12 

96. 33 (+167) 

3.58 (-46) 

13 

89. 92 (+75) 

8.  84 (+147) 

14 

94.47 (-18) 

5.  27 (-56) 

15 

93. 30(-103) 

6. 1 3 (+1 72 ) 

16 

87.51 (+171) 

10. 8 3 (— 22 ) 

17 

88. 53 (+101) 

9. 68 (-150) 

18 

92. 30(+63) 

6. 54 (+153) 

19 

92.08 (+172) 

7. 36(-30) 

20 

87 . 04 (+76) 

11. 22 (+147) 

21 

88. 57 (-25) 

10. 22 (-59) 

22 

86.40(-122) 

12. 09(+121) 

23 

86. 02 (+160) 

12. 02 (-56) 

24 

90. 69 (+69) 

8.  56 (+1 39 ) 

25 

90. 50 (+27) 

8.45 (+36) 

26 

89.11(+4) 

8. 96(+18) 

27 

79. 53(+104) 

16. 56 (-167) 

28 

85. 7  9 (+2 ) 

12. 17  (-24 ) 

29 

75. 67 (+88) 

18. 92 (+140) 

30 

82. 85(+159) 

14.46(— 72) 

31 

79. 36(+71) 

17. 1 5  (+1 26  ) 

32 

82. 51(+10) 

14 . 65 (+12) 

33 

91. 50 (-17) 

7 . 87 (-46) 

34 

88. 99 (—17) 

9. 51 (-55) 

Second 
Harmonic 
0.25  sec 

Third 

Harmonic 

0. 19  sec 

Fourth 
Harmonic 
0.15  sec 

0. 27 (-107 ) 

0. 66  (-137) 

0. 00 (-34 ) 

0.  98(+89) 

0. 06(-107) 

0. 07 (-127 ) 

0. 1 3 (+153) 

0.  01 (-1 75) 

0. 10(-13) 

0. 03 (+114) 

0. 02 (-134) 
0. 03 (+1 11) 
0.00(-100) 
0. 03 (-126) 
0.00(+30) 

1. 24 (-105) 

0. 41 (-51) 

0. 31 (+172) 
0.51(-74) 

0. 76 (+122) 

0. 21 (-21) 

0. 05(-117) 

0. 03 (+129) 

0. 07 (+1 55) 

0. 09 (-80) 

0.07 (+48) 
0. 02 (-141 
0. 00(-l 25) 
0. 00(+15) 

0. 01 (+105) 

1. 10(+90) 
0.07 (+103) 
1.03(-136) 

0. 22 (-100) 
0.51(+101) 

0. 22 (-6) 

0. 02(+81) 

0. 1 6 (— 62 ) 
0.04 (+124) 

0.  06 (+36) 

0. 08 (-115) 
0.01 (-100) 
0. 05 (-7) 
0.0K+103) 
0.00(-100) 

1. 40(+150) 

1 . 47 (-36) 

0. 97  (-115) 

0. 52 (+126) 
1.45  (— 14  5 ) 

0. 17 (-48) 

0.  27 (+75) 

0.  17 (-28) 

0.  03 (-102) 

0.  20(-80) 

0. 08(+l 1 ) 

0. 05(+159) 
0. 02 (+65) 
0.02 (+42) 
0. 09(-24) 

1 . 06 (-88) 

1. 38 (+4) 
1.57 (+88) 

0. 68 (-138) 

0. 95 (+52) 

0. 11(-124) 

0. 11(-120) 

0. 33(+115) 

0. 04 (-48) 

0. 09(+72) 

0. 04 (-165) 
0.02 (+88) 
0. 06(+32) 
0. 02 (-21 ) 
0.02 (+56) 

1.61 (+14) 

3. 11(-81) 
1.68 (-51) 

4. 05(+19) 

2. 19 (+54) 

0. 29 (+18) 
0.57  (-4 ) 

0.  30 (— 1 06) 

0.  95 (-114) 
0.41 (+160) 

0. 03(+14) 
0.24 (+66) 

0. 06(-139) 
0.41 (+105) 
0. 10(-90) 

2. 85(+176) 

2.  2 9 (+1 9 ) 

0.  61 (- 72 ) 

1 . 32 (-91 ) 

0.48(-144) 

0.43(+25) 

0.  02(-65) 

0. 17(-117) 

0. 18 (-114 ) 
0.13 (+18) 

0. 00(+114) 
0.01 (-166) 

*  Ocean  gage. 


Table  11 


Distribution  of  Wave  Energy  Throughout  the  Model,  20-deg  Incident 


Direction,  Period  =  1.00  sec.  Ocean  Wave  Height  =*  0.112  ft 


Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Gage 

Fundamental 

Period 

1.00  sec 

First 

Harmonic 

0. 50  sec 

Second 

Harmonic 

0- 33  sec 

Third 
Harmonic 
0.25  sec 

Fourth 
Harmonic 
0.20  sec 

1* 

99 . 06 (-177) 

0 . 74 (+57) 

0 . 20(-104) 

0 . 00 (+71) 

0.00(+57) 

2* 

98. 31 (-154) 

1 . 58(+78) 

0.11(-148) 

0.01(+86) 

0.00(-25) 

3 

98. 57 (+170) 

1. 38(-70) 

0.04(+67) 

0 . 00 (-123) 

0.00(-95) 

4 

99. 49 (+109) 

0. 50(-88) 

0.01(+159) 

0.00 (+24) 

0. 00(-122) 

5 

97.40(+55) 

2.57 (+116) 

0.03(+130) 

0.00 (+67) 

0.00(+141) 

6 

98.95(+2) 

0. 98 (—33 ) 

0 . 07 (-72) 

0.00(-73) 

0.00(-120) 

7 

99. 77 (-57) 

0 . 22 (-94) 

0.00(-107) 

0.00 (+164) 

0 . 00 (+1 31 ) 

8 

99.08(-112) 

0. 84(+154) 

0 . 07(+88) 

0.01 (+25) 

0.00(-30) 

9 

97. 80(-158) 

2. 15(+47) 

0.03(+113) 

0 . 01 (-76) 

0.00(+118) 

10 

99.41 (-172) 

0. 54(+19) 

0.05(-65) 

0.01 (+146) 

0.00(-20) 

11 

95. 33(-16) 

4 . 54 (-29) 

0. 12  (-22) 

0 . 01 (+45) 

0.00(+47) 

12 

97. 14 (—65 ) 

2 . 82 (-165) 

0. 04(+136) 

0.00(-142) 

0.00(+178) 

13 

99. 53(-122) 

0.44 (+122) 

0 . 02 (+113) 

O.Ol(-O) 

0.00(-117) 

14 

97.05(-177) 

2. 88 (+5) 

0.06(-123) 

0.01 (+130) 

0. 00(-68) 

15 

96.99(+122) 

2. 99(-131) 

0 . 01 (-68) 

0 . 00 (-46) 

0.00(+61) 

16 

98. 88(+71) 

1.07 (+107) 

0.04(+121) 

0.00(+151) 

0.00(+178) 

17 

99 . 55 (+34) 

0. 31 (+79 ) 

0. 13(-140) 

0.01(-112) 

0 . 00  (-140) 

18 

99. 33(+21) 

0. 63 (+42) 

0.04(+115) 

0. 00(+172) 

0 . 00 (+4 3) 

19 

97 . 58(-179) 

2. 40(+6) 

0 . 02 (+163) 

O.OK-95) 

0.00(+128) 

20 

92. 19 (+131) 

7.47 (-112) 

0. 33(+27) 

0. 01 (+130) 

0.00(-41) 

21 

94 . 65 (+76) 

5 . 24 (+120) 

0. 11 (+179) 

0 . 01 (-166) 

0 . 00 (-5) 

22 

96 . 85 (+13) 

3.09 (+15) 

0.04 (-18) 

0.01(-61) 

0.00(-51) 

23 

96. 87 (-47) 

2.91 (-90) 

0. 20(-133) 

0 . 02 (+114) 

0.00(+82) 

24 

95. 20(-97) 

4 . 76 (+160) 

0.04(+86) 

0. 00 (+95) 

0. 00(+7) 

25 

97. 88(-127) 

2.09(+85) 

0.03(-25) 

0.00(-79) 

0.00(-128) 

26 

98. 14(-132) 

1.81 (+27) 

0.05(+114) 

0.00(-74) 

0.00(+23) 

27 

86. 53 (+40) 

11.47 (+49) 

1 . 67(+69) 

0. 28 (+9  5 ) 

0 .06(+l 35) 

28 

87.48(-12) 

11. 14(-41) 

1 . 21 (-66) 

0. 15 (-83) 

0. 02(-76) 

29 

82. 13(-136) 

15.57 (+120) 

1 . 95 (+114) 

0.31 (+42) 

x - ^ 

vC 

r—l 

1 

■<T 

c 

c 

30 

84 . 25 (-1 32) 

13. 39 (+70) 

1 . 94 (-80) 

0. 36(+131) 

0.06(+4) 

31 

91 . 61 (+171) 

7 . 65 (-49) 

0.67 (+111 ) 

0.06(-82) 

0.02 (+107) 

32 

94. 48 (+12  7 ) 

5. 14 (-118) 

0. 36  (-1 2) 

0.01 (+89) 

0.00(-l 52) 

33 

94. 7  3 (+103) 

5 . 18(— 170) 

0. 09(-50) 

0.00 (-154) 

o 

x— ' V 

1 

■O 

> — r 

34 

96. 77(+98) 

3 . 19 (+1 79) 

0. 01 (+79) 

0.03 (+1 71) 

0.00(-l 20) 

*  Ocean  gage. 


Table  12 


Distribution  of  Wave  Energy  Throughout  the Model,  20-deg  Incident 
Direction,  Period  =  1.00  sec,  Ocean  Wave  Height  =  0.134  ft 
Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Fundamental 

First 

Second 

Th  i  rd 

Fourth 

Period 

Harmonic 

Harmonic 

Harmonic 

Harmon i c 

Cage 

1.00  sec 

0.50  sec 

0.33  sec 

0.25  sec 

0.20  sec 

1* 

96.87 (+11 ) 

3. 09 (+32) 

0.01 (-124) 

0.03(-167) 

0.01 (+177) 

2* 

97 . 32 (+36) 

2. 43(+99) 

0. 20 (+9 9) 

0.03(+151) 

0.02 (-131) 

3 

98 . 84 (-8) 

1.14 (+13) 

0 . 00 (-5) 

0.01 (-144) 

0 . 00(+165) 

4 

95. ll(-67) 

4 . 69 (- 1 37) 

0.17 (+1 34) 

0.02 (+68) 

0.01 (+39) 

5 

96.96(-121) 

2. 73 (+81) 

0. 28(-75) 

0.03(+144) 

0.00 (+41) 

6 

99 . 61 (-17  3) 

0. 33 (-22) 

0.04 (-168) 

0 . 02 (+15) 

0.00 (+12 5) 

7 

98.01 (+129) 

1 . 82 (-91) 

0 . 14 (+102) 

0. 02 (-89) 

0.01(+b8) 

8 

96. 51(+74) 

3. 32 (+160) 

0. 14 (-145) 

0.02 (-86) 

0.00 (+4 7) 

9 

96.17^+31) 

3 . 82 (+38) 

0.00 (+104) 

0.00(-80) 

0.00 (-6) 

10 

98.94 (+19) 

1.0 3 (+9) 

0.03(+34) 

0.00 (+150) 

0.00 (-108) 

11 

97.01 (+169) 

2.81 (-54) 

0. 1 5 (+1 14 ) 

0.02 (-50) 

0.00 (+160) 

12 

97 . 44 (+114) 

2 . 42 (-123) 

0.11(+7) 

0 . 02 (+145) 

0.00 (-7 7) 

13 

93.71 (+59) 

5. 82(+129) 

0.41 (-152) 

0 . 04 (-63) 

0 . 02 (+29) 

14 

91 . 82 (+7 ) 

7. 76(+l) 

0. 38(-5) 

0.03(-12) 

0.01 (+40) 

13 

96. 50(-54) 

3. 16 (-126) 

0 . 29 (+158) 

0.05 (+94) 

0 . 00(+66) 

16 

98. 63(-104) 

1. 33C+130) 

0 . 02 (-136) 

0.01 (+106) 

0.00 (+4) 

17 

98.51(-139) 

1. 37 (+105) 

0. 1 0 (— 3  5 ) 

0.02 (-177) 

0 .01 (+9 1 ) 

18 

98.00(-148) 

1.95 (+67) 

0.05 (-130) 

0.00(+71) 

0. 00(-81 ) 

19 

87.69 (+3) 

11.28 (+0) 

0. 90(+6) 

0.11 (+31 ) 

0.0 3 (+76) 

20 

90 . 00 (-48) 

9 . 08 (-127) 

0 . 84 (+174) 

0.08(+121) 

0.01 (+91) 

21 

94 . 78  (-107) 

5.06(+127) 

0. l6(+0) 

0.00(-101) 

0.00(-79) 

22 

91. 33(-167) 

7. 90 (+20) 

0. 70 (— 1 53) 

0.07 (+41) 

O.Ol(-lll) 

23 

89 . 84 (+135) 

9.01 (-97) 

0 . 9  9  (  +4  2  ) 

0. 14 (-1 70) 

0.02(-4) 

24 

91.14 (+87) 

7 . 79 (+153) 

0. 9  3(-l 35) 

0. 1 3 (—53) 

0.02 (+63) 

25 

96 . 92 (+60) 

2 . 96 (+96) 

0 . 10 (+85) 

0.02 (+97) 

0.00 (+159) 

26 

96. 79 (+58) 

3.05(+92) 

0.1 5 (+80) 

0.01 (+66) 

0.00 (+14 7) 

27 

82. 30(-144) 

14 . 74 (+58) 

2 . 51 (-90) 

0.40(+128) 

0.06(+3) 

28 

74 . 67 (+165) 

19.91(-54) 

4.46(+99) 

0. 84 (- 103) 

0.1  3 (+56) 

29 

74.71 (-88) 

19 . 99 (+117) 

4.41 (-88) 

0.80 (+4) 

0. 09(+88) 

30 

77.16 (+43) 

1 7 . 85 (+60) 

4 . 04 (+83) 

0. 79 (+109) 

0 . 1  6  (  + 1  2  5  ) 

31 

84.00(-14) 

13.54 (-47) 

2 . 14 (-66) 

0. 28(-77) 

0.03(-6l) 

32 

84. 02 (-54) 

13. 33(— 125) 

2 . 19 (+1 79) 

0. 38 (+128) 

0.08 (+92) 

33 

88. 28(-72) 

10. 27 (—166) 

1  . 24 ( +1 14) 

0.17 (+39) 

0.04 (-20) 

34 

91 . 48(-74) 

8. 09 (-163) 

0 . 40(+l 15) 

0.02 (+48) 

0.01 (-22) 

*  Ocean  page. 


Table  13 

Distribution  of  Wave  Energy  Throughout  the  Model,  20-deg  Incident 


Direction,  Period  =  1.50  sec,  Ocean  Wave  Height  -  0.077  ft 
Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Ca8e 

Fundamental 

Period 

1.50  sec 

First 
Harmonic 
0.75  sec 

Second 
Harmonic 
0.50  sec 

Third 
Harmonic 
0.38  sec 

Fourth 

Harmonic 

0 . 30  sec 

1* 

98. 68 (+45) 

0. 97 (-132) 

0 . 15  (-38) 

0.19(+34) 

0.00(-10) 

2* 

98. 85(+59) 

1.01  (-69) 

0.09 (+4 7) 

0.01 (+157) 

0.04 (-72) 

3 

91 . 28 (+175) 

8. 34(-19) 

0. 3  3 (+128) 

0. 05(-l 79) 

0.00(-58) 

4 

90. 32 (+142) 

9. 27 (-98) 

0. 38(-36) 

0.03(-84) 

0.00(-113) 

5 

92.05(+108) 

7 .  78  (-167) 

0. 1 2 (— 1 02 ) 

0.05(-10) 

0. 00(+176) 

6 

9 3. 4 3 (+81) 

6. 47(+118) 

0.06(-130) 

0.02 (+134) 

0.01(-153) 

7 

93. 79 (+48) 

6.00(+42) 

0. 14(+117) 

0.04 (-111) 

0.03(-17) 

8 

95. 70 (+15) 

3. 79 (-32) 

0 . 32(+10) 

0 . 16 (-18) 

0.03(-169) 

9 

98. 82 (-4) 

1 . 10(-86) 

0.07 (+154) 

0.0K-112) 

0 ,00(-l 79) 

10 

99. 18  (-1 1 ) 

0. 80 (-79) 

0.00(+19) 

0.01 (+45) 

0.00(+163) 

11 

91. 70(-71) 

7 . 88(-149) 

0.40(-160) 

0.01  (-27) 

0 .02 (-1 26) 

12 

90 . 74 (-105) 

8. 28(+143) 

0 . 09 (+72) 

0.04(+110) 

0 . 00(-35) 

13 

90.93(-139) 

8.40(+72) 

0.  56 (—53) 

0. 12 (+170) 

0.00(-165) 

14 

90. 22(-170) 

9. 13 (+4) 

0. 55(-178) 

0.09 (+10) 

0.02(-166) 

15 

89. 52(+158) 

9 . 51 (-70) 

0. 82(+72) 

0.07 (+1 36) 

0.08(-33) 

16 

90. 07 (+127) 

9 . 82 (-145) 

0.08(-75) 

0 . 02 (+111) 

0.01 (+170) 

17 

94 . 89 (+110) 

4 . 73 (+180) 

0 . 37 (-47) 

0.01 (-6) 

0.00 (+44) 

18 

96.08(+112) 

3. 86 (-x78) 

0 . 04 (-82) 

0.01 (+102) 

0.01 (+73) 

19 

94. 16(+54) 

5. 64(+115) 

0.11(+144) 

0 . 05 (-87) 

0.03(— 1 11 ) 

20 

90.67 (+23) 

9 . 07 (+47 ) 

0. 21 (+62) 

0. 04 (+23) 

0.02 (-32) 

21 

91 . 95 (-16) 

7. 69 (-34) 

0. 35(-92) 

0. 00(-5) 

0.00 (+160) 

22 

88. 48(-48) 

10.91(-110) 

0.56(-170) 

0.05(-176) 

0 . 00 (+161 ) 

23 

86. 68(-82) 

12. 50(-179) 

0.57 (+125) 

0. 16(-0) 

0.09 (-50) 

24 

88. 54 (-107) 

10. 32 (+128) 

1.06(+1) 

0.57 (-62) 

0.29(+l 38) 

25 

92.99(-120) 

6.95(+93) 

0 . 02 (-32) 

0 . 03 (+41) 

0 .01 (-56) 

26 

94. 38 (-119) 

5 . 49 (+97) 

0.11  (-6) 

0.02 (-169) 

0.0] (+10) 

27 

93.45(-165) 

5 . 82 (+8) 

0. 57(+178) 

0. 1 5 (+78) 

0.0] (-72) 

28 

90 . 99 (+1 59) 

7. 84 (-63) 

0.88(+79) 

0.27 (-126) 

0.03(-6) 

29 

88. 58(+120) 

8.  54(-142) 

2 . 6  6  ( -  3  4 ) 

0. l8(+74) 

0.04(+l 57) 

30 

86.21 (+87) 

11 .06(+160) 

2.  SO ( —  1 35) 

0. 19 (-62) 

0.04(+l 7) 

31 

85. 17 (+57 ) 

13. 36(+102) 

1 . 2  8 (+ 1 2 1 ) 

0.  14 (  +  1 62) 

0 . 06 (- 1 28 ) 

32 

85 . 93(+38) 

12. 33 (+S1 ) 

1 .47(+88) 

0. 2(5  (+12  7) 

0.01 (+104) 

33 

90. 25 (+29) 

8. 97(+22) 

0. 69(+29) 

0.09(+78) 

0.00 (-38) 

34 

90 . 83 (+2  7 ) 

8.61 (+21) 

0. 5  3 (+21 ) 

0. 03 (- 3) 

0.00 (+26) 

*  Ocean  gage. 


Table  14 

Distribution  of  Wave  Energy  Throughout  the  Model,  20-deg  Incident 
Direction,  Period  =  1.50  sec,  Ocean  Wave  Height  =  0.111  ft 
Percent  of  Total  Energy  and  Phase  Angle  (deg) 


Gage 

Fundamental 

Period 

1.50  sec 

First 
Harmonic 
0.75  sec 

1* 

96.99(-154) 

2. 52(-161) 

2* 

97. 55(-139) 

1.93(-119) 

3 

84.22(-22) 

14.41(-67) 

4 

85.16(-57) 

14.42(-153) 

5 

86. 80(-87) 

12. 68(+138) 

6 

87 . 21  (-119) 

11. 87 (+61) 

7 

87.95(-152) 

10. 68(-5) 

8 

94.15(+177) 

5. 38(-83) 

9 

98. 79(+156) 

1 . 00 (-137) 

10 

98. 90 (+150) 

1.00(-132) 

11 

82. 66 (+86) 

16. 15(+167) 

12 

76.94(+54) 

20. 20(+98) 

13 

79.42(+20) 

17.81(+26) 

14 

77 . 50 (-12) 

20. 01 (-49) 

15 

77.99 (-43) 

20 . 46 (-121) 

16 

88. 53 (-68) 

10.25(+172) 

17 

93.01(-83) 

6. 42 (+137) 

18 

93. 5  7 (—87 ) 

6. 13(+133) 

19 

86. 74(-154) 

11. 64 (+63) 

20 

79. 36(+175) 

17. 80(-16) 

21 

80. 85 (+138) 

15 . 84(-96) 

22 

73. 87 (+107) 

20. 35 (-163) 

23 

72. 27 (+79) 

22.49(+133) 

24 

81 . 15 (+56) 

17. 46(+70) 

25 

85. 88 (+44) 

12. 70(+39) 

26 

83. 7  5 (+44 ) 

14.41(+41) 

27 

84. 51 (-13) 

12. 06(-50) 

28 

77.52 (-48) 

16.14 (-118) 

29 

74. 19(-88) 

19 . 04 (+166) 

30 

66. 54 (-121) 

23. 29 (+101) 

31 

65. 34 (-145) 

25.17 (+44) 

32 

69. 88(-166) 

20. 88(-0) 

33 

71.98(— 175) 

23.47 (-27) 

34 

74.98(-175) 

21 . 05 (-22) 

Second 
Harmonic 
0.50  sec 

Third 
Harmonic 
0.38  sec 

Fourth 
Harmonic 
0. 30  sec 

0.40(+55) 
0.47 (+105) 

1 . 30(-115) 

0. 30(+125) 

0. 50(+47) 

0.06(-65) 
0.05 (+2) 
0.06 (+124) 
0.08 (+40) 
0.01(+138) 

0.03(+177) 
0.00(+151) 
0.0K+118) 
0.03(-26) 
0.00 (+6 7) 

0.7 3 (-60) 
1.27(-169) 
0.46(+59) 
0.12 (+120) 
0.05(+4) 

0. 16(+176) 
0.07 (-12) 
0.00(-105) 

0. 07 (-21) 
0.05(-176) 

0.02 (+49) 

0.03(+161) 

0.00(+122) 

0.01(+127) 

0.00(+ll) 

1 .12  (-89) 

2. 36(+143) 
2.36(+30) 

2. 29 (—87 ) 

1. 36(— 175) 

0.05 (-55) 

0. 32 (-159) 

0. 33(+23) 

0. 19(-96) 

0. 16(+119) 

0.0K+45) 
0.09(-105) 
0 . 08(+47) 
0.02(-135) 
0.03(+72) 

1. 12(+108) 

0. 52(+13) 

0. 29 (+1) 

1 . 42(-99) 

2. 35(+147) 

0.09 (+40) 
0.04(-94) 

0. 00 (-131) 
0.17 (+88) 
0.42 (-27) 

0 . 00 (-21) 

0 . 01 (+138) 
0.00 (+2 9) 
0.03(-53) 

0 . 08 (+153) 

2. 67 (+40) 

4 . 75 (-69) 

4 . 11 (-162) 

1 . 36(+101) 
1.29 (+83) 

0 . 54 (+169) 

0. 86(+37) 
0.92(-98) 

0 . 28 (+130) 
0.12 (+109) 

0.09(-50) 

0 . 1 8 (+1 34 ) 
0.21(-25) 
0.0K+112) 
0.0K+124) 

1 . 72 (+74) 

2 . 72 (-70) 

4 . 69 (+172) 
4.77 (+62) 

7. 59 (-40) 

0.11(+75) 

0. 54(-97) 

1 . 30(+110) 

1 . 54 (-49) 
2.05(-172) 

0.01 (+90) 

0. 17(-119) 
0. 36 (+5 3) 
0.46(-147) 
0. 53 (+59 ) 

6. 80(-118) 

7 . 49 (+166) 

3. 57 (+1 31) 

3. 31 (+138) 

2.04 (+82) 

1 . 35 (-18) 

0. 78 (-65) 

0. 58(-57) 

0. 65(-73) 

0. 39 (+161) 
0 . 19 (+94) 
0.08(+110) 

*  Ocean  gage. 
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Photo  4.  Typical  wave  pattern  for  1,5-sec  wave  of  height  0.083  ft  near  the  wave  generator 
approaching  from  an  incident  direction  of  20  deg  (upwave  and  dovmwave  region) 


Photo  5.  Typical  wave  pattern  for  1.5-sec  wave  of  height  0,079  ft  near  the  wave  generator 
approaching  from  an  incident  direction  of  20  deg  (dovnwave  region) 


approaching  from  an  incident  direction  of  20  deg  (upwave  and  downwave  region) 


mpw  n’n'i  wit 


for  1, 0-sec  wave  of  heigh 
ident  direction  of  20  deg 


Photo  11.  Typical  wave  pattern  for  0.75-sec  wave  of  height  0.176  ft  near  the  wave  generator 
approaching  from  an  incident  direction  of  20  deg  (upwave  and  downwave  region) 


Photo  12.  Typical  wave  pattern  for  0.75-sec  wave  of  height  0.176 


Photo  15.  Typical  wave  pattern  for  0.75-sec  wave  of  height  0.123  ft  near  the  wave  generator 
approaching  from  an  incident  direction  of  20  deg  (downwave  region) 


movement,  time  after  starting  generator  =  0.5  min,  wave  period  *  1.0  sec 
height  near  generator  *  0.139  ft,  angle  of  incidence  *  30  deg 


:lme  after  starting  generator 
generator  *  0,139  ft,  angle  < 


Photo  21.  Met  dye  novenent,  tine  after  starting  generator  ■  4.0  min,  wave  period  “  1.0  sec 
wave  height  near  generator  ■  0.139  ft.  angle  of  Incidence  *  30  deg 


DIFFRPCT ION  THEORY  DISTPNCE  PLONG  BREPKWPTER  =12.00  FT 


PLATE  7 
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PLATE  10 


PLATE  12 


DIFFRACTION  THEORY  DISTANCE  ALONG  BREAKWATER  =6.00  FT 


PLATE  14 


I FFRACT I  ON  THEORY  DISTANCE  ALONG  BREAKWATER  =8.00  FT 


PLATE  15 


DIFFRfiCTION  THEORY  DJSTPNCE  PLONG  BREPKWPTER  =10.00  FT 


PLATE  18 


DISTANCE  ALONG  BREAKWATER  =6.00  FT 


DISTANCE  ALONG 


Ho  =  0. 122  FT 

DISTfiNCE  PLONG  BREPKWPTER  =12.00  FT 


Ho  =  0. 180  FT 

DISTANCE  ALONG  BREAKWATER  =8.00  FT 
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PLATE  28 
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PLATE  31 
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PLATE  34 
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D1STRNCE  PL0NG  BREPKWPTER  =8.00  FT 


PLATE  35 


DISTANCE  ALONG  BREAKWATER  =10.00  FT 


DISTRNCE  PLONG  BREPKWPTER  =10.00  FT 


DISTANCE  ALONG  BREAKWATER  =12.00  FT 


Ho  =  0.077  FT 

DISTPNCE  PLONG  BREPKWPTER  =6.00  FT 


DISTANCE  ALONG  BREAKWATER 


DISTANCE  ALONG  BREAKWATER  = 8 . OG  FT 


PLATE  47 


DISTANCE  ALONG  BREAKWATER 
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PLATE  51 


AVERAGE  BOTTOM  VELOCITIES 

- >  COMBINED  REFRACTION  AND  DIFFRACTION 

FT/SEC  ANGLE  OF  INCIDENCE  =  30  DEGREES 

WAVE  PERIOD  =  1.0  SECONDS 
WAVE  HEIGHT  NEAR  GENERATOR  =  0.139  FEET 


APPENDIX  A:  NOTATION 


bI,1>2 


Wave  amplitude,  ft 

Second  derivative  of  wave  amplitude  in  x-direction,  1/ft 

Second  derivative  of  wave  amplitude  in  y-direction,  1/ft 

Wave  ray  spacing,  ft 

Wave  speed  of  propagation,  ft/sec 

Wave  group  velocity,  ft/sec 

Friction  coefficient,  dimensionless 

Local  water  depth,  ft 

Water  depth  at  breaking,  ft 

Local  rate  of  energy  dissipation,  lb/ft/sec 

Base  of  system  of  natural  logarithms,  2.71828,  dimensionless 

2 

Wave  energy  density  at  breaking,  ft-lb/ft 

2 

Wave  energy  density  in  deep  water,  ft-lb/ft 

Mathematical  function  for  wave  equation  in  cylindrical 
coordinates ,  dimensionless 

Wave  thrust,  lb/ft 

2 

Gravitational  constant,  ft/sec 

Local  wave  height,  ft 

Wave  height  at  breaking,  ft 

Incident  wave  height,  ft 
1/2 

(-1)  ,  dimensionless 

Wave  number,  2tt/L,  1/ft 
Diffraction  coefficient,  dimensionless 
Refraction  coefficient,  dimensionless 

Combined  refraction  and  diffraction  coefficient,  dimensionless 
Wavelength,  ft 

Ratio  of  group  velocity  to  wave  velocity,  dimensionless 
Longshore  current  parameter,  dimensionless 
Rate  of  wave  power  transmission,  Ib/sec 

Longshore  current  parameter  from  Longuet-Higgins,  dimensionless 
Radius  in  cylindrical  coordinates,  ft 
Bottom  beach  slope,  dimensionless 


♦ 


A1 


T  Wave  period,  sec 

u  Instantaneous  particle  velocity 

u  Two-dimensional  vector 

u  Parameter  combining  breaker  wave  energy  density  and  local  water 
depth,  ft/sec 

u  Maximum  wave  orbital  velocity,  ft/sec 

ul,u2  Diffraction  integration  limits,  dimensionless 

V^g  Longshore  current  velocity,  ft /sec 

x  Distance  offshore,  ft 

a  Local  wave  angle  with  shoreline  or  breakwater,  deg 

a.  Breaker  angle  of  waves  with  shoreline,  deg 

D 

Deepwater  wave  approach  angle,  deg 

6  Angle  between  wave  ray  and  breakwater,  deg 

2 

6  Defined  parameter,  (a  +  a  )/k  a,  dimensionless 
r  xx  yy 

Ab  Wave  ray  spacing,  ft 

ri  Local  water-surface  elevation,  ft;  complex  wave  amplitude 

0  Direction  in  cylindrical  coordinates,  deg 

0  Angle  between  breakwater  and  undiffracted  wave  crests  beyond 
breakwater,  deg 

2 

y  Horizontal  eddy  viscosity  coefficient,  lb-sec/ft 

e  2  4 

p  Fluid  density,  lb-sec  /ft 

2 

t  Local  shear  stress,  lb/ft 

2 

T-  Local  mean  longshore  shear  stress,  lb/ft 
Is  2 

<p  Velocity  potential,  ft  /sec 

u)  Wave  angular  frequency,  2tt/T,  l/sec 

8  Partial  derivative  symbol,  dimensionless 

V  Horizontal  gradient  operator,  dimensionless 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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